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The NLO Brodsky-Lepage-Mackenzie (BLM) scale for the pion transition form factor has been 
determined. To achieve that, a consistent calculation up to n/-proportional NNLO contributions to 
both the hard-scattering amplitude and the perturbatively calculable part of the pion distribution 
amplitude has been performed. By combining and matching the results obtained for these two 
amplitudes, a proper cancellation of coUinear singularities has been established and the 75 ambiguity 
problem (related to the use of the dimensional regularization method) has been resolved by using 
the naive-75 as well as the 't Hooft-Veltman (HV) schemes. It has been demonstrated that the 
prediction for the pion transition form factor is independent of the factorization scale /^l^ at every 
order in the strong coupling constant, making it possible to use the simplest choice /if? = Q^ at the 
intermediate steps of the calculation. Assuming the pion asymptotic distribution amplitude and 
working in the MS scheme, we have found the BLM scale to be ^\ — fj.%]^nj ~ Q^/9. Based on 
the same distribution, the complete NLO prediction for the pion transition form factor has been 
calculated in the av definition of the QCD coupling renormalized at fj,'^ = fiy 
It is in good agreement with the presently available experimental data. 

PACS numbers: 13.40.Gp,12.38.Bx,13.60.Le,11.15.Me,13.40.Hq 
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I. INTRODUCTION 

The pion transition form factor plays a crucial role 
in testing QCD predictions for exclusive processes. It 
appears in the amplitude that relates two, in general 
virtual, photons with the lightest hadron, the pion, 
7*7* -^ tt". Historically, this process attracted much 
interest since the axial anomaly ||l[ fixes the value of the 
form factor when both virtualities of the photons are zero 
(77 ^ tt"). For large virtualities of the photons (or at 
least for one of them) perturbative QCD (pQCD) is ap- 
plicable g] . A specific feature of this process is that the 
leading-order (LO) prediction is zeroth order in the QCD 
coupling constant, and one expects that pQCD for this 
process may work at accessible values of spacelike pho- 
ton virtualities H. Experimentally, the most favorable 
situation is when one of the photons is real (7*7 — > tt''). 

The framework for analyzing exclusive processes at 
large-momentum transfer within the context of pQCD 
was initiated and developed in the late seventies P, El. 
It was demonstrated to all orders in perturbation thc- 



Alexander von Humboldt Fellow. On leave of absence from the 
RiiHjer Rosknvir TTis titiitp Zagreb, Croatia.; Electronic address: 



melic@thphys.irb. hi 



T Electronic address: iizic@thDhvs.irb.hr 



-^Electronic address: aassekOthphys. irb.hr 



ory that exclusive amplitudes at large-momentum trans- 
fer factorize into a convolution of a process-dependent 
and perturbatively calculable hard-scattering amplitude, 
with a process- independent distribution amplitude (DA) , 
one for each hadron involved in the amplitude. Whereas 
the DA is intrinsically nonperturbative and its form is 
determined by some nonperturbative methods, the DA 
evolution is subject to a perturbative treatment. In the 
standard hard-scattering approach (sHSA), a hadron is 
regarded as consisting only of valence Fock states, trans- 
verse quark momenta are neglected (collinear approxima- 
tion) as well as quark masses. 

Although the LO predictions in the sHSA (as well as in 
the modified hard-scattering approach (niHSA) in which 
the collinear approximation is abandoned B) have been 
obtained for many exclusive processes, only a few pro- 
cesses have been analyzed at the next-to-leading order 
(NLO): the pion electromagnetic form factor 0, ^, |j, 
the pion transition form factor |7|^ (and [^ in the 
mHSA), and the process 77 -^ MM {M = tt, K) ^. 

It is well known that, unlike in QED, one cannot rely 
upon the LO predictions in pQCD (the expansion param- 
eter, i.e., the running coupling constant is rather large at 
current energies), and that higher-order corrections are 
important. The size of the NLO correction as well as 
the size of the expansion parameter, i.e. the QCD run- 
ning coupling constant, can serve as sensible indicators 
of the convergence of the expansion. However, as the 



truncation of the perturbative expansion at any finite or- 
der causes the residual dependence of the prediction on 
the choice of the renormahzation scale and scheme, these 
choices introduce an ambiguity in the interpretation of 
the finite-order perturbative prediction. In general, in- 
cluding higher-order corrections has a stabilizing effect 
(see H , for illustration) reducing the dependence of the 
predictions on the schemes and scales (since the all-order 
prediction is independent of the scheme and scale choice) . 
However, to asses the convergence of the perturbative ex- 
pansion, it is necessary not only to extend the calculation 
beyond the LO (which is a very demanding task in many 
cases), but also to optimize the choices of the scale and 
scheme according to some sensible criteria. 

In the Brodsky-Lepage-Mackenzie (BLM) procedure 
[ [13| , |lj], all vacuum-polarization effects from the QCD 
/3-function are resummed into the running coupling con- 
stant. Since the coefficients /3o,/3i,--- are functions of 
n/ (number of flavors), according to BLM procedure, the 
renormahzation scale best suited to a particular process 
in a given order can be, in practice, determined by com- 
puting vacuum-polarization insertions in the diagrams 
of that order, and by setting the scale demanding that 
rij-proportional terms should vanish. The renormahza- 
tion scales in the BLM method are physical in the sense 
that they reflect the mean virtuality of the gluon prop- 
agators and the important advantage of this method is 
"pre-summing" the large {Poas)"" terms, i.e., the infrared 
renormalons associated with the coupling constant renor- 
mahzation ( |15|| and references therein). 

The optimization of the renormahzation scale and 
scheme for exclusive processes by employing the BLM 
scale fixing was elaborated in Ref. |1^ . It was stated that 
exclusive processes are especially sensitive to the choice 
of the renormahzation scale for the underlying hard- 
scattering amplitude and since each external momentum 
entering an exclusive reaction is partitioned among many 
propagators of the underlying hard-scattering amplitude, 
the physical scales that control these processes are in- 
evitably much softer than the overall momentum trans- 
fer. The BLM method was applied to the pion electro- 
magnetic form factor and the 77 — > tt+tt^ process. For 
the pion transition form factor, the size of the BLM scale 
was only assumed (taken the same as for the pion electro- 
magnetic form factor). Since the LO prediction for the 
pion transition form factor is zeroth order in the QCD 
couphng constant, the NLO corrections |0, g, 10 rep- 
resent leading QCD corrections and the vacuum polar- 
ization contributions appearing at the next-to-next-to- 
leading order (NNLO) are necessary for determining the 
BLM scale for this process. 

The purpose of this work is to determine the BLM scale 
for the pion transition form factor, i.e., for the 7*7 — > tt 
process. Although, the structure of the process is simple, 
the calculation of higher-order corrections to the hard- 
scattering amplitude is complicated by the 75 ambigu- 
ity, which appears when using dimensional regulariza- 
tion. In our calculation we use dimensional regulariza- 



tion in D = 4 — 2e dimensions to regularize both ultravi- 
olet (UV) and coUinear singularities. We have obtained 
the LO, NLO, and n/-proportional NNLO terms for the 
hard-scattering amplitude using the Feynman gauge and 
modified minimal-subtraction scheme (MS) for which a 
suitable compact form has been adopted. In order to 
correctly subtract the coUinear singularities and also to 
verify the right choice of the 75 prescription, we have also 
determined the LO, NLO, and n/-proportional NNLO 
terms of the perturbatively calculable (evolutional) part 
of the distribution amplitude. The 75 ambiguity present 
in the calculation of the hard-scattering part has been 
resolved by combining and matching the results for the 
hard-scattering amplitude with the results for the distri- 
bution amplitude part. The proper cancellation of sin- 
gularities has been established and the 75 problem has 
been resolved using both the so-called naive-75 [|l6[ and 
the 't Hooft-Veltman (HV) schemes 0, |l8) . It has been 
demonstrated that the prediction for the pion transition 
form factor is independent of the factorization scale ^'^ at 
every order in strong coupling constant. Finally, we have 
been able to justify the natural choice fij^ = Q^ for the 
factorization scale and to determine the renormahzation 
scale using the BLM scale setting method. 

The plan of the paper is as follows. Section II is de- 
voted to some preliminary considerations. In Sec. Ill the 
calculational procedure is briefly outlined. In Sec. IV the 
LO, NLO, and the ^/-proportional NNLO unrenormal- 
izcd contribution to both the hard-scattering amplitude 
and the perturbatively calculable part of the pion distri- 
bution amplitude are obtained. Renormahzation of the 
UV divergences and factorization of the coUinear diver- 
gences present in the hard-scattering and the pion distri- 
bution amplitude are performed in Sec. V. The results 
for both amplitudes are obtained in the naive-75 as well 
as the HV scheme. The complete leading-twist analytical 
expression for the pion transition form factor up to Uf- 
proportional NNLO terms is obtained in Sec. VI. Section 
VII is devoted to determining the BLM scale for the pion 
transition form factor based on which the complete NLO 
numerical predictions are then obtained in the AIS and 
ay renormahzation schemes. The concluding remarks 
are given in Sec. VII. The 75 problem is addressed in 
detail in Appendix A. In Appendix B the Feynman rules 
for the perturbatively calculable part of the distribution 
amplitude are derived. Finally, in Appendix C, we clar- 
ify some often obscured points on the coupling constant 
renormahzation and justify our renormahzation conven- 
tion. 



II. PRELIMINARIES 



The pion transition form factor F-y*-y7r(Q^) for a pseu- 
doscalar meson tt" is defined in terms of the amplitude 

r^" for 7*(g, ^) -I- 7(fc, v) -^ 7r(P), as 

F'^" =ze^ F^.^^(Q^) £'^""'5 P^qp . (2.1) 
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FIG. 1: Feynman diagram describing the 7*7 ^ qq amplitude 
in terms of which the hard-scattering ampHtude for the 7*7 —> 
IT transition is obtained. 



For large- momentum transfer Q^ ~ —q^ , the form factor 
can be represented [pi W as a convolution 



F^,^^{Q^) = ^*{x,^il) 



Th{x,Q^,h],), 



(2.2) 



where 
by 



stands for the usual convolution symbol defined 



(2.3) 



A{z)(^B{z) = / dzA{z)B{z) 



In Eq. (2^), the function Th {x , Q"^ , fij^) is the hard- 
scattering amplitude for producing a coUinear qq pair 
from the initial photon pair; $*(x,/i|,) is the pion dis- 
tribution amplitude representing the amplitude for the 
final state qq to fuse into a pion, i.e., the probability am- 
plitude for finding the valence qq Fock state in the final 
pion with the constituents carrying fractions x and (1 — 2:) 
of the meson's total momentum P; /x|, is the factoriza- 
tion (or separation) scale at which soft and hard physics 
factorize. 

The hard-scattering amplitude Th is obtained by eval- 
uating the 7*7 — ^ qq amplitude, which is described by 
the Feynman diagrams in Fig. |l], with massless on-shell 
quarks coUinear with outgoing meson. By definition, T^ 
is free of coUinear singularities and has a well-defined 
expansion in asipji), with /i|j being the renormalization 
(or coupling constant) scale of the hard-scattering ampli- 
tude. Thus, one can write 



TH{x,Q'^,^ip) 

.(0) 



T^^>^,,Q^) + ^^Tl,'\x,Q^f,l) 



«s(Mfi) ^(2) 



(47r)2 



n^'ix,Q',fij„fi%) + --- . (2.4) 



evolution kernel 
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^^Viix,u) 
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V2{X,U) 



(2.6) 



If the distribution amplitude $(x, /i§) is determined at an 
initial momentum scale /ig (using some nonperturbative 
meth ods) , then the differential- integral evolution equa- 
tion ( |2.5| ) can be integrated using the moment method 
to give ^{x^^'^p) at any momentum scale /x|n > /Zq. The 
one-Q and two- loop |1^, ^, |l| corrections to the evo- 
lution kernel are known, but because of the complicated 
structure of the two- loop corrections, it was possible to 
obtain numerically only the first few moments of the evo- 
lution kernel [E2l. However, based on the conformal spin 
expansion, the conformal Ward identities, and the con- 
formal consistency relation, the complete analyti cal form 
of the NLO solution of the evolution equation ( |2.5| ) has 
been obtained in Rcf. [ p3| . 

Instead of using ^{x^^Pp), one often introduces the dis- 
tribution amplitude (/)(x,/i|n) normalized to unity 



dx(t){x, fip) = 1 



and related to ^{x,iJ,jp) by 

^ix,fil) = Ni 
where 

N^ = 



{x,^J.f) 



fn 



(2.7) 



(2., 



(2.9) 



is the normalization constant imposed by the leptonic 
7r+ -^ fJ-^i^fi decay, /^r = 0.131 GeV is the pion decay 
constant, and Nc{— 3) is the number of QCD colors. 

It is convenient to expand the distribution amplitude 
0(x,/Zq) (determined at the initial scale /Iq) in terms of 
the Gegenbauer polynomials Cj {2x — 1), representing 
the eigenfunctions of the LO evolution kernel Vi : 



cj){x,fil) =6xil-x) J2'BnCf/\2x--l). (2.10) 



n=0 



The nonperturbative input is now contained in the -B„ 
coefficients and ^ ' denotes the sum over even indices. 
The DA (2.10), when evoluted to the scale ^|,, is repre- 
sented by the perturbative expansion 



(a;,AiF) = '/' (a^,MF) + 



asip-p) 



An 



cb''-^ix,f,'p) + ^ 



(2.11) 



where 



,LO 



Although the function $(a;, ^^) is intrinsically nonper- 
turbative (containing the effects of confinement, nonper- 
turbative interactions, and meson bound-state dynam- 
ics) , it satisfies an evolution equation of the form 

d 

^|^-^$(x,/i|.) = y(x,M,as(Ai|^)) (g) $(u,^|,), (2.5) ^NLO 

where V{x,u,as{fJ''^)) is the perturbatively calculable 



{x,fip) = 6a::(l - x) 



{x,fj.p) = 6x(l - x) 



Y,'Bl:^{^j?p)Cl'\2x~l). 

n=0 
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\t^ 



)C'J\2x 



(2.12) 

!)■ 
(2.13) 



The coefficients B^^^{fij,) and B^^^{n%) depend on the 
nonperturbative input Bn, as well as on the scales /Iq 
and fjijp. Their exact form can be read from the results 
obtained from |23l and listed in [R|^. 



III. CALCULATIONAL PROCEDURE 

Before proceeding with the calculation, we would like 
to point out some subtleties connected with the calcula- 
tional procedure that (we think) deserve more explana- 
tion. 

The hard-scattering amplitude Th is obtained by eval- 
uating the 7* + 7 -^ qq amplitude, which contains 
coUinear singularities, owing to the fact that final state 
quarks are taken to be massless and on-shell. Since, by 
definition, Th is a finite quantity free of coUinear sin- 
gularities, these singularities should be subtracted. The 
7* -|- 7 — > gg amplitude with the Lorentz structure fac- 
tored out as in Eq. (2.1) and denoted by T factorizes 



T{u,Q^) ^THix,Q'^,^i%) (g) ZT,coiix,u;n%), (3.1) 

where, as usual, u and I — u denote the quark/ antiquark 
longitudinal momentum fractions, /i|^ is a factorization 
scale at which the separation of coUinear singularities 
takes place and all coUinear singularities are factorized 
in Zt^coI- 

On the other hand, a process-independent distribution 
amplitude for a pion in a frame where P+ = P'^ + P'^ = 1, 
P- =P° - P^ = 0, and Pi = is defined E, Eol, M as 



$(u) 



where 



dz 
X ( 



^j(„_(l_„))j;-/2 






TT ) , (3.2) 

(Z+=Z_L=0) 



fl = exp<ig / dsA^{zs)z /2 > 



(3.3) 



is a path-ordered factor making $ gauge invariant. The 
matrix element in (3^) contains an ultraviolet diver- 
gence coming from the light-cone singularity at z^ = 
Ig, Q. This divergence should be regulated, and af- 
ter renormalization, which introduces a renormalization 
scale /i|., z^ is effectively smeared over a region of or- 
der z^ — — z^ ^ ^/i^'b.- As 9- result, the pion distri- 
bution amplitude $(u,/i|j) is obtained corresponding to 
the pion wave function integrated over the pion intrinsic 



transverse momentum up to the scale jjj^. The distribu- 
tion amplitude <i>(w,/i^) is a fi nite quantity and enters 
the convolution expression (2.2). 

The un reno rmalized pion distribution amplitude <i>(u) 
given in (3.2) and the distribution amplitude ^(w,/!^) 
renormalized at the scale /2|j are (owning to the mul- 
tiplicative renormalizability of the composite operator 
^7+75$!^) related by a multiplicative renormalizability 
equation 



$(u) = Z^,reniu,V;flj^)(g)^{v,flj^) 



(3.4) 



By differentiating this equation with respect to p,^ one 
obtains the evolution equation ( |2.5| ), with the evolution 
potential V given by 



V 



^Z2 



r-,2 



d 



(3.5) 



For notational simplicity, here and where appropriate, we 
use the notation in which the convolution ((g)) is replaced 
by the matrix multiplication in x-y space (unit matrix is 

given by 1 = S{x — y))- 

The pion distribution amplitude as given in ( |3.2| ) , with 
Jtt) being the physical pion state, of course, cannot be de- 
termined using perturbation theory. If the meson state 
Jtt) is replaced by a \qq;t) state composed of a free 
(coUinear, massless, and on-shell) quark and antiquark 
(carrying momenta tP and (1 — i)P), the amplitude (3.2) 
becomes 



(t>iu,t) 



dz 
X (0 



J{u-{l-u))z- /I 



vI/(-z)l-Jn*(z) 



qq;t 



(3.6) 



Taking (3.6) into account, Eq. ( |3.2| ) can be written in 
the form 



$(u) =: (j){u,t)® {qq;t\Tr) 



(3.7) 



The amplitude (3.6) can be treated perturbatively, mak- 
ing it possible to investigate the high-energy tail of the 
pion DA, to obtain Z^^ren and to determine the DA evo- 
lution. 

The distribution 0(m, t) is multiplicatively renormaliz- 
able and the UV singularities that are not removed by the 
renormalization of the fields and by the coupling constant 
renormalization, factorize in the renormalization con- 
stant Ziji^ren at the renormalization scale jl ^. A part from 
UV singularities, the matrix element in ( |3.6[) contains 
coUinear singularities (since the initial state quarks are, 
as before, taken to be massless and on-shell), which are 
absorbed in Z^,.oi at the factorization scale /Zq. Hence, 
one obtains 



^ It should, however, be pointed out th at, in contrast to the ex- 
pansion parameter a5(/i|,)/(47r) from (2.1]), the expansion pa- 
rameter a.s{fi^p)/Tr was chosen in |q]. Hence, the expressions for 

Bn from PI should be modified accordingly. 



(u,i) 

= Z(li.ren{u,V]jl\)> 



(jjviv, s; fiji, /ig) ® Z^,coi (s, t; fi^) . 

(3.8) 
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FIG. 2: Pictorial representation of the of the pion transition 
form factor calculational ingredients: T represents the pertur- 
batively calculable hard-scattering am plitu de , while $ is the 
pion distribution amplitude given by (3.2) which can be ex- 
pressed, as in (3.7), in terms of the perturbatively calculable 
part (j) (3^) and the perturbatively uncalculable part. 



the distribution 



Upon combining Eqs. (3.7) and (3 
<^{u) can be written in the form 

$(u) = Z^^ren{u,V;flj^) (g) 0y (u, S; /i|j, ^q) ® $(s, /ig) . 

(3.9) 



Here, 



<^{s,fil) = Z^^coi{s,t;nl)(E) {qq;t\TT) 



(3.10) 



represents the nonperturbative input (containing all ef- 
fects of coUinear singularities, confinement, and pion 
bound-state dynamics) determined at the scale /Uq, while 
4>v(y, s; fijj^, /Iq) governs the evolution of distribution am- 
plitude to the scale jij^: 

"^{v, lA) = <^v(w, s; A?j, Mo) ® ^{s, ^2) ^ 



(3.11) 



and satisfies the evolution equation 



r-.2 



d 



P-R^7r^4)v{v,s,jj.ji,HQ) 



•d^i 



Viv,s',^l) ® (l)v{s',s,p.%til). (3.12) 



By convoluting the "unrcnormalized" (in the sense 
of coUinear singularities) hard-scattering amplitude 
T{u, Q^) with the unrcnorm alize d pio n di stribution am- 
plitude <&(m), given by Eqs. ( |3.lD and (3A), respectively, 
one obtains (in a way analogous to |g, |25|] ) the following 
expression for the pion transition form factor F^t^T^iQ'^)- 

F^,^^{Q^) = $t(^) ^ T{u, Q2) . (3.13) 

The divergences of T(u, Q^) and <i>(u) cancel for flj^ = ii\ 

Zt,coi{x, u; ^i%) (g) Z^^reniu, v] ^%) = S{x - v) , (3.14) 



and the usual expression (|2.2D emerges, where the pion 
transition form factor is expressed in a form of the con- 
volution of two finite amplitudes: 



It is worth pointing out that the scale /i|. representin g the 
boundary between the low- and high-energy parts in ( |2.2| ) 
is, at same time, the separation scale for coUinear singu- 
larities in T(u, Q^), on the one hand, and the renormal- 
ization scale for UV singularities appearing in the per- 
turbatively calculable part of the distribution amplitude 
<&(«), on the other hand. The calculational procedure 
explained above is illustrated in Fig. |[ 

Our main goal in this work is to determine the BLM 
scale for the pion transition form factor. To achieve 
that, we make use of the calculational procedure out- 
lined above and in the following sections calculate the 
LO, NLO, and n/-proportional NNLO contributions to 
the perturbative expansions of both the hard-scattering 
amplitude and the distribution amplitude. 

This is the first calculation of the hard-scattering am- 
plitude T{u, Q^) of an exclusive process with the NNLO 
terms taken into account. The subtraction (separation) 
of coUinear divergences at the NNLO is significantly more 
demanding than that at the NLO. Owing to the fact 
that the process under consideration contains one pseu- 
doscalar meson, the calculation is further complicated 
by the 75 ambiguity related to the use of the dinrensional 
regularization method to treat divergences. 

In order to correctly subtract the coUinear divergences 
and determine the right treatment of the 75 matrix, we 
determine the LO, NLO, and n/-proportional NNLO 
contributions to the distribution (f>{u,t) given in (p^), 
and by that, following (3.7-3.9|), the renormalization con- 



stant Z^^ren and the evolutional part (jyy of the distribu- 



tion amplitude <I>(w) (3.2). Since, there is no 75 ambigu- 
ity in the DA calculation, the 75 ambiguity present in the 
hard-scattering calculation is resolved using (3.14). As an 
additional check, we employ two 75 schemes in our cal- 
culation. Finally, we obtain the (75 scheme independent) 
prediction for the pion transition form factor F^*^.,^ up 
to the n/-proportional NNLO terms, expressed in terms 
of the finite quantities Th{x,Q^, fij^) and ^{x,fip). 



IV. LO, NLO, AND ^/-PROPORTIONAL NNLO 

UNRENORMALIZED CONTRIBUTIONS TO 

THE HARD-SCATTERING AND THE 

DISTRIBUTION AMPLITUDES 

In this section we present the calculation of the LO, 
NLO, and n/-proportional NNLO contributions to the 
hard-scattering amplitude and the perturbatively calcu- 
lable part of the distribution amplitude. 



A. Contributions to the hard-scattering amplitude 

The hard-scattering amplitude Th for the pion transi- 
tion form factor is obtained by evaluating the 7*7 — > qq 
anrplitudc for the parton subprocess, which is described 
by the Feynman diagrams of Fig. nl 

The qq pair has to be projected into a negative-parity 
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FIG. 3: Lowest-order Feynman diagrams contributing to the 
7*7 ^ W amplitude. 
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and spin (pseudoscalar) state. This is achieved by in- 
troducing the projection operator 75 /P/\/2 and taking 
the trace over a fermion loop. On the other hand, the 
color-singlet nature of the qq state is taken into account 
by introducing the factor X]q=i ^ap/VN^, and taking the 
trace over the color indices. Also, the flavor function 
{uu~ dd)/V^ should be included. 

The hard-scattering amplitude T{u, Q^) can generally 
be expressed as an expansion in as 






(47r)2 
+ •••} 



(4.1) 



LO contributions 



In the LO approximation there are only two Feynman 
diagrams that contribute to the 7*7 — > qq amplitude. 
They are shown in Fig.^. The contribution of diagram A 
(after ie^e'^''°'f^Paqp is factored out) is given by 



Nt 



1 



0^ 1 



where 



and 



Nt 



Gtt 



^/2 



V2 
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(4.2) 



(4.3) 



(4.4) 



is the factor taking into account the flavor content of the 
qq pair. The contribution of diagram B is obtained by 
making the replacement u — > (1 — u) in ( [4.2[ ). Therefore, 
the lowest-order (QED) contribution to the 7*7 -^ qq 
amplitude , i.e., to T{u,Q'^) given in (4T), is 



1 



1 — M 



(4.5) 



FIG. 4: Distinct one-loop Feynman diagrams contributing to 
the 7*7 -^ qq amplitude, generated from diagram A of Fig. H 
by inserting a gluon line. 



2. NLO contributions 

At NLO there are 12 one-loop Feynman diagrams con- 
tributing to the 7*7 — > qq amplitude. They can be gen- 
erated by inserting an internal gluon line into the lowest- 
order diagrams of Fig. |^. We use the notation where 
Aij is the diagram obtained from diagram A by insert- 
ing the gluon line connecting the lines i and j, where 
i,j — 1,2,3. Since NLO, and all higher-order, diagrams 
generated from diagram B can be obtained from the cor- 
responding diagrams generated from diagram A by using 
the substitution u ^ (1 — w), the total number of NLO 
diagrams to be evaluated is 6. They are shown in Fig. |j. 

These diagrams contain ultraviolet (UV) singularities, 
and owing to the fact that the final state quarks are mass- 
less and on shell they also contain coUinear singularities. 
To regularize these singularities, we use dimensional reg- 
ularization in _D = 4 — 2e space-time dimensions. 

As is well known, dimensional regularization leads to 
an ambiguity when dealing with the pseudoscalar matrix 
75. The reason for this lies in the fact that the matrix 75 
cannot be unambiguously defined in D ^ 4 dimensions. 
In practice, the difficulty arises in evaluating a trace con- 
taining a single 75. We address this problem in detail in 
Appendix ^. 

In order to make sure that our results for the pion 
transition form factor are 75 scheme independent, we 
have evaluated all the contributions using two schemes: 
the naive-75 scheme p6| and the 't Hooft-Veltman (HV) 
scheme [|l^ |l^, defined by (Al) and (A3), respectively. 

A few remarks concerning the diagrams with quark 
self-energy corrections where the quark momentum p is 
on-shell {All, A22) are in order. Since these corrections 
modify external legs, each of these diagrams is accompa- 
nied with a factor of 1/2 coming from the expansion of 
the quark field renormalization constant ^/Z2. In dimen- 
sional regularization, the contributions of each of these 
diagrams turn out to be proportional to (p^) and, 
therefore, vanish when p^ = 0. On closer inspection. 



TABLE I: Contributions T^^ and T^^ (defined in ([4.q)) of Feynman diagrams sliown in Fig. W. Apart from tlie result denoted 
by 125=1, the listed results correspond to the results obtained using the HV scheme, or equivalently, to the results obtained in 
the naive-75 scheme with 75 being outside the contracting 7 matrices. The result denoted by 12^=1 has been obtained using 
the naive-75 scheme with 75 placed between the contracting 7 matrices and between two Dirac slashed loop momenta. 



«J 



T 



Aij 



l-Aij 



11(22) 

33 

23 

13 
12a=o 
125=1 



-i(l-)(l-|)(l-^)- 



-(l-6)(l-n)- 

(1-.)(1-^)-^ 

(1-.) 





u 







-i(l-)(l-|)(l-^)- 



(2 + .)(l-n)-^ 



[2 + ,){l~u)-^+{^^+e^Ul-{l-uD 

-2(1 + .) Q + 2) l—!i{l-[l-uD 

_2(l-e)ii^i^(l-(l-u)-^) 
e u 



however, one finds that this vanishing is a result of the 
cancellation of a UV pole with coUinear pole. The UV 
pole contributes to the renornialization of the quark fields 
(already taken into account by the factor V^) and even- 
tually leads to a correct running of the coupling constant. 
The contribution of any of the diagrams Aij shown in 
Fig. m can be generally expressed as 



'^Aij — -7TT 



Nt 



1 



1 



u An ■' 



(4.6a) 



where Cf = 4/3 is the color factor (the same for all 
diagrams), while Taij is defined by 



TAij — 



.(0) 
UV 



(e) 



1 



l-2e 



pUV 
'-Aij 



+r™(.)r^r« g, , («b) 



with the following abbreviations 

^uvi^) - Me) r(l-2e) 



{47ry , (4.7a) 



rfi(.) =r(l + .)MXl_l)(4,).. (4.7b) 



The first F function on the right-hand side of Eqs. ( |4.7| ) 
originates from the loop momentum integration, while 
the integration over Feynman parameters produces Fs 
collected in a fraction. Conse quen tly, the singularity con- 
tained in F(e) appearing in (4.7a) is of UV origin, while 
the singularity contained in F(— e) appearing in ( 4.7h| ) is 
of infrared (IR) origin. It should be pointed out, how- 
ever, that none of the diagrams of Fig. contains a soft 
(genuine IR) singularity, so that the, here and in the fol- 
lowing, subscript (and / or the superscript) IR signifies 



the coUinear singularity. If the relation 



F(z)F(l-z) = - 



is taken into account in (4.7), one finds that 



.(0) 



(0), 



T'uUe) = -rY>ie) 



(4.8) 



(4.9) 



Nevertheless, we continue to keep track of the origin of 
the UV and coUinear singularities. 

The contributions T^^ and T^^- of the individual di- 
agrams are given in Table |l Following the explanations 
and notation given in App. |A|, we list the contributions 
obtained using the HV scheme, which are equivalent to 
the results obtained in the naive-75 scheme with 75 be- 
ing positioned outside the contracted 7 matrices. For 
diagram A12 we also list the contribution obtained in 
the naive-75 scheme corresponding to the case where 7 
matrices are contracted through the string of 7 matrices 
of the form /75 /. As elaborated in Appendix H, the 75 
ambiguity in diagrams All, A22, A33, A23, and A13 has 
been resolved with the help of QED Ward identities. The 
remaining ambiguity in diagram Al2 is parameterized by 
6, taking the value for the first choice for handling 75 in 
diagram A12, and 1 for the second. Our results listed in 
Table | are in agreement with jl^ (but see the comments 
in App. 0). 

The NLO contribution r(i)(-u) from Eq. (Q is of the 
form 



r(i)(u) = 



r[?Ue)A:TW(-) 



l-2e 



.(0) 



(1), 



+rrii(e)Y3^nii(«) 



^ ) ^ (4-lOa) 



v/\yVv\A>- 






"^^ 








^^^ 




^ 






1 


^AAyW\i * 


J . 



FIG. 5: Distinct vacuum-polarization two-loop Feynman di- 
agrams contributing to the 7*7 — > qq amplitude, which have 
been obtained from the one-loop diagrams of Fig. H by in- 
serting the vacuum-polarization bubbles. 



where 



7p(l) / X 

TuvW 
= Cf 



1 



1 — M 

+ {u -^ 1 - u) , 



|(i-^r + i(i-(i-.r) 



(4.10b) 



and 



7p(l)/ N 



= Cp 



1-u 



4 + 2e-[l^le+^e^]{l-uy 



i'2 , „ , 1 — u , 
+ - - 4 - e + 4,5 - 3 2 -I- e 



(l-(l -«)"'=) +{u^l-u). (4.10c) 



5. Uf -proportional NNLO contributions 

By inserting the vacuum polarization bubbles in the 
NLO diagrams of Fig. H, the NNLO diagrams displayed 
in Fig. are obtained. The vacuum polarization inser- 
tion is (in the Feynman gauge) given by the replacement 



-*5 



where 



kX ^ab 



P + irj 



-« .9 



5ab 



P + irj J P + irj 



(4.11) 



Il{p) = 



(Att) 



IT] 



1 



(5-3e)--(l-e)n; 

-r(^)'^Vr^'l\"^^4.) 



(l-|e)(l-2e) " r(l-2e) 

(4.12) 
and, due to gauge invariance of the complete finite order 
contribution (which we have used as an additional check 
of our calculation), it can effectively be described by 



~*5 



6ab 



P + irj 



-*.9 



kA "ab 



P + irj 



n{P 



(4.13) 



We are interested only in n/-proportional part (from the 
quark loops inserted in the gluon propagator) 



n„,(P) 



where fnj is defined by 



fnf{e,fi^) = 



[P + ir]Y (47r)- 



: fnf{e,fi^), 



(4.14a) 



1 -e 



l-irfj (l-|e)(l-2e) 



T^k^ 



(4.14b) 



The contributions of the two-loop Feynman diagrams 
shown in Fig. a can than be generally written as 



HAtj)r 



where 



Nt 1 






T, 



(4.15a) 



T, 






(1-6) 



^(0) (^ p(l) / X [^ ~gj ^UV 

,2\ 2e 



while, similarly to (4.7), the abbreviations 



.(1) 
uv 



(e) = 



r(2e) r(l-2e)r(l-e) 



r(l + e) r(l-3e) 



(4.15b) 



(47r)'^ (4.16£ 



have been introduced. 

The contributions T, ..'.-, of the individual diagrams 
are listed in Table |l[ The integrals Ii{u,e) {i = 1,2,3) 
appearing in this table are defined as 



/i(u, e) = I{u;e,2e) 
h{u,e) = /(w;e,l + 2e) 
h{u,t) = /(u; 1 + e, 1 -f 2e) , 



(4.17) 
(4.18) 
(4.19) 



where 



I{u;a,c) = / dy- 

Jo (1 - uy) 



1 + a 



2i^i(c,l-Ka,2 + a;u). (4.20) 



As far as the 75-scheme dependence of the NNLO dia- 
gram contributions is concerned, it is the same as for the 
NLO diagrams from Table B. 



TABLE 11: Same as Table fn but for contributions TFAij) f ^^^ (Aij) f (defined in (WTSl)), corresponding to the Feynman 
diagrams shown in Fig. fel 



y 






rpiR 



l (l-e)-(f-f)' 
2(l + f)(l-H 



11(22) 



33 



23 



13 



12^=0 



1 (1-6)^(1-1) 

2(l + f)(l-|.) 



(1-u)- 



(1-e 



(l-|e) 
{1-ef 



3-v(l --)"'' 



(1-1^) 




3 ^(1-^)"'^ 



3 ^ h{u,e) 



{i-uY 



(2 - e - 2^2 



(1-1^) 



(1-^)- 



(1-1^) 



[26(1 - 6) hiu, 6) - (2 - 6 - 26^) U J3(U, e)] 



-2(1 + 6) [-(1 - u)-''" - 2eh{u, 6) - (1 - e)uh{u, e)] (1 - u) 



12^=1 







-2(1 - 6) [-(1 - u)-^' - (1 - e)uh{u, 6)] (1 - u) 



(2)/ 



The n/-proportional NNLO contribution Tn/ (u) from Eq. ( [I.l| ) takes the form 



where 



T(2."/)(m) = 



(l-|e)(l-26)(l-3e) 



r[;^(e)r^^|.(e), 



Tuv (") 



-(0) ^.M-(l)/ 



(1-^) 



+rr/y(e)r,;,(e)^^_2^^^^_2e)(l_3e) 



7p(2,«/)/ X 



Q2 



2£ 



77^(2,"/) 



Tuv m = C'i= 



1 (1-^)^ 

l-u 1 - ie 



-•- + 2 



— (2,n 



T^/fl, (w) = C'F- 



l+ ie- Se"- T^e 



9^2 _ 13 
2" 2^ 2 



-2e 



(l-e)2 



1-if 



(1 



(l + i)(l-|e) 
■-{2S-2){l + e){l-u) 



+ 2(1- (2(5-l)e)(l-u) 
/2(u,e) 



-f (u ^ 1 — m) , 



1-i^ 



- l-(2<5-l)6 (l-e)(l-y) 



+ {u — > 1 — m) . 



uh{u,e) 



(4.21a) 



(4.21b) 



(4.21c) 



Our results are expressed in a compact form in which 
the complete functional dependence on the dimensional 
parameter e is retained. This is in contrast to the expan- 
sion over e often encountered in the literature. In this 
expansion, nonleading terms in e are neglected before 
renormalization and factorization of collinear singulari- 
ties. As we show in Sec. M, it is advantageous (both for 
the simplicity and accuracy check of the calculation) not 
to expand the functions ^[/yjji over e. 



B. Contributions to the perturbatively calculable 
part of the distribution amplitude 

111 Se c. II we have defined the distribution amplitude 
4>{u,t) (3^) representing the perturbatively calculable 
part of the pion distribution amplitude. Following po| , 
we have rederived the Feynman rules for this operator in 
the Feynman gauge. They are listed in Appendix ^ 

We now proceed to calculate LO, NLO, and Uf- 
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tp 




(l-t)p 



FIG. 6; The LO diagram contributing to the 4>{u, t) distribu- 
tion (O), i.e., (Q. 






proportional NNLO contri buti ons to the 0(u, t) distribu- 
tion amplitude defined in (|3.6|), or equivalently in (B3). 



Contrary to Efl, M, we use dimensional regulariza- 
tion to regularize both UV and coUinear singularities ^ . 
This enables us to combine these results with the hard- 
scattering results, also obtained by employing the dimen- 
sional regularization. Compared to the hard-scattering 
amplitude calculation, calculation of the (j) amplitude is 
complicated by the f act that nonc ovariant Z"*", l~ and 
(5- function terms (see (B4) and (B5)) appear in the loop- 
momenta. To deal with these types of terms and in or- 
der to simphfy the expressions we follow the prescription 
given in |20| . 

The presence of two 75 matrices in the traces has 
enabled their unambiguous treatment in the naive-75 
scheme (see Appendix ^). Additionally, we have ob- 
tained the results using the HV scheme, which, however, 
introduces the "spurious" anomalous terms, and hence 
the additional rcnormalization is required. The corre- 
sponding rcnormalization constant will be determined by 
comparing the results obtained in the naive-75 with those 
obtained in the HV scheme. 



The perturbatively calculable i 
represented as a series in as 



{u, t) amplitude can be 



{u,t) = 0(o)(M,i)-|-^0(i'(u,i) 



(4^)2 



47r 



(4.22) 



^ The evolutional behavior of the DA can be extracted from (3.6) 
even when using dimensional regularization for both UV and 
mass singularities. We introduce the auxiliary scale Q^ and we 
insist on discriminating between UV and collincar singularities. 
Otherwise, the UV and coUinear part of higher-order corrections 
would cancel, leading to <f>{u, t) = S{u — t). 





FIG. 7: The one-loop diagrams contributing to the (j>{u, t) 
distribution (3.(3), i.e., (B3) (in the Feynman gauge). 



1. LO contribution 

The contribution of the LO Feynman diagram dis- 
played in Fig. ^ reads 



a(o) 



(M,t) ^5[u~t). 



(4.23) 



2. NLO contributions 



At NLO there are 5 one-loop Feynman diagrams con- 
tributing when the Feynman gauge is used. They are 
displayed in Fig. 0. The general form of these individual 
contributions (denoted by X) is given by 



(px{u,t) 

__ as 
47r 

where 



Kx{u,t) 



(4^)2 



d^l 



1 



(27r)^ [P + i-qf 



Kx{u,t)^2Cf 



dyFxiu,t;y), 



(4.24) 



(4.25) 



and the Fxju, t;y) and Kx{u,t) contributions can be 
read from Table III. The notation ui = u, W2 == 1 — u 



has been introduced , as well as the usual definition of 
the "+" form 

[f{x,v)] =F{x,y)~5{x-y) I dzF[z,y), (4.26) 
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TABLE III: Contributions Fx appearing in Eqs. (4.25) and (4.31) and correspo nding to th e dia grams {X) displayed in Fig. W 
and Fig. H. We list also the general results for the integrals appearing in Eqs. (4.25) and (4.31), and parameterized by k = 
and K = e, respectively. 



Fx{u,t;y) 



dyy'^Fx{u,t;y) 



A 



A + BI + B2 



CI 



(1 - t) [u\5{ui_ -yti) + u2S{u2 - y i2)] 



(1 - e) I ui(5(mi —yti)+ U25{u2 -yt2)\ 



Ul 



iti - ti 



5{ui - yti) 



(1-e) 






a-.H|.„.-«, I 



+ 



Itl -^ U2 
tl -^t2 

Ml -^ U2 
tl -^t2 



tl ii-^1 ' \ti 



the presence of which reflects the chiral symmetry con- 
servation. This "-I-" form is a consequence of the fact 
that the axial current is conserved in the chiral limit, 
and represents a general all-order property pl[|. 



By definition, the D dimensional integral in Eq. (4.24) 
gives zero in dimensional regularization, but only if we do 
not distinguish between UV and coUinear singularities '^. 
By discriminating between the singularities of different 



^ The D dimensional integrals appearing in Eqs. (4.24) and (1.30) 
are of the form 

d^l 1 



(27r)-D («2 +ir;)° 



By employin] 



d^l 



{l-pf 



(27r)-D (^2 + iri)°'{(l - p)2 + iri) 



VO 



and insisting on distinguishing the F functions obtained from the 
loop-momentum integrations and V functions from the Feynman 
parameter integration, it can be shown that 



1 



(47r)2 (p2 + ir?)«- 
"F(a-2-|-e) F(3 



4-71/^2 



■ a - 



e)V{2 - e) 



+ 



F(a) F(5-Q:-2e) 

F(a - 1 -I- e) F(2 - « - e)r(3 - 



(2-.) 



r(a) 



F(5- 



■2e) 



Here, the first fraction in the terms containing F functions cor- 
responds to the loop-momentum integration possibly resulting 
in UV singularities, while the second fraction corresponds to 
the integration over Feynman parameters and consequently, to 
coUinear singularities. For a < 2 + e, where e = or <g, only the 
UV singularities appear, while for a > 2 -\- e, only the coUinear 
singularities appear. The two terms in the bracket cancel in both 



cases, so I\ 



[a#2 + 6],f« 



: 0. However, for a = 2 + e, both UV and 



coUinear singularities are present, and the cancellation can occur 
only if we abandon distinguishing them. 



origin, we obtain the following expression 

d^l 1 " 



-i (47r) 



Ai 



{2tt)D {p + zr,)^ 



UV 



(e) 



1 



l-2e 



'^-(^^T^O-f) 






(4.27) 



where Q- > represents the auxiliary scale. 

The NLO contribution (j)^^^{u,t), to which Feynman 
diagrams of Fig. [7| contribute, can then be expressed by 



(!)(., i) = 



.(0) 
UV 



(^) 



l-2e V 2/ ^ ' ' 



+^S<^'Ti 0-1) ''■'"<-«] (I 

(4.28) 

where the function K^^' calculated in the naive-75 scheme 
amounts to 

1 



kW 




e{t - u) 



(4.29) 



3. Uf -proportional NNLO contributions 

By insertingthe vacuum polarizations in the NLO dia- 
grams of Fig. Yi we obtain the NNLO diagrams displayed 
in Fig. H. The n/-proportional contributions'* of these 
diagrams, evaluated using Eqs. ( 4.13 - 4.14 ), read 



* There are more two-loop diagrams containing quark loops, but 
they contribute to the ny-proportional NNLO part of the meson 
singlet distribution amplitude. 
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FIG. 8: Distinct two-loop Fey nman diag rams contributing 
to the 4>{u,t) distribution (3.6), i.e., (B3), which have been 
obtained by inserting the vacuum polarization bubbles in the 
one-loop diagrams of Fig. M. 



The function /„j,(e, ^^) is defined in ( 4.14b| ), while 
K(x)„f is given by 



K^x)^^{u,t)^2CF / dyy'Fxiu,t;y), (4.31) 

Ja 



and Fx{u,t;y) and Ki^x),^j{u,t) can be read from Table 

ill. 



Si milarl y to (4.27), the D dimensional integral from 
Eq. (|4.30|) gives 



-i {Anf 



d^l 



1 



2 \ ' 






l-e _ 1 



^^^^'\l-3e)(l-|6)l + 



+r\%) 



l-e 



1 



(l-3e)(l- §e) 1 + e 



(4.32) 



(f>(x)„f{u,t) 



(47r)2 



(-^^v) ^(X)„,(^,i)(l + e) ,fn,(e,M') 



(47r) 



d^l 



{2tt)D (P + iT]) 



2+e 



(4.30) 



where Q^ > represents the auxiliary scale. 

The n/-proportional NNLO contribution of the dia- 
grams displayed in Fig. M takes the form 



4>(^'-f\u,t) 



r(o) felr^'^ (e) 



(1-6)2(1-1) 



3.,a-2e)(l-3e)(l-|e) 

, r(o) (e)T^^he) (i-£)'(i~i) 



(1-|6)(1-26)(1-36)(1-|6) 
I 



X(2,«/)(u,t) 



^ 

Q' 



(4.33) 



where the function K'-^''^f^ calculated in the naive-75 
scheme amounts to 



K'-^'''f\u,t) = 2Cf 



t 



(l-e) + 



t — u 



u — > 1 — u 



e{t-u) 



(4.34) 



4- The HV scheme results 

The preceding results have been calculated in the 
naive-75 scheme. When the HV scheme is used, only the 
results for A diagram from Fig. |7| and the corresponding 
"bubble" diagram from Fig. S differ from the naive-75 
results, and are given by 



= 4>A{u,t) (l + Ae 



(l-f)(l-e) 



, (4.35) 
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and 



<^(A)„/"'0 = '/'(A)„/(u,0 + ^'/'(A)„/ (■«,*) 



respectively. These results, obviously, these results bear 
the signature of chiral symmetry violation. 



0(^)„,Kt)(^l + 4e^^-— ^-^ 



(4.36) 

Hence, when using the HV sc heme , the functions K'-^^ 
and /i:(2,«y) j^ gqg ([4.28|) and (t4.33|) get replaced by 



'(1) _ 



K)^>, = 2CFUe 



1 



-9{t~^u) 



u ^ 1 — u 



+ K'^^\ (4.37) 



and 



K 



(2,"/) 
HV 



{u,t) 



2 Cf <^ 4e 



1-f 
u ^ 1 — u 
f ^ 1 -t 



t(t)'«"- 



V. RENORMALIZATION AND 

FACTORIZATION OF COLLINEAR 

SINGULARITIES 



A. General renormalization procedure 



Since in this work we present the calculation (up to Uf- 
proportional NNLO contributions) of the hard-scattering 
amplitude T{u, Q^) (4.1), as well as of the perturbatively 



calculable part of the DA 4){u,t) (4.22), both containing 
, UV and coUinear singularities, here we outline the gen- 

eral procedure for the renormalization of UV and the 
(4.38) factorization of coUinear singularities. 

I 



We introduce the amplitude M (having the same form of the perturbative expansion as the amplitudes T(u, Q^) 
and (j){u,t)): 



47r (47r)2 



(-lnf)M^^'-f^ + --- 



(5.1a) 



where 



A^W - 



rZ{e){a^^ + ear + e'ar+0{e^)) 



Tf^ie){4'^ + eai^ + e'4^ + 0{e')) 



EL 



A^(^-/) = [r[,^(e)r«(.)(C^^%.6;' 



n„ay_^^2^rv,C/V^Q^^3) 



, p(0) 



(e)r^^i(e) (6[,'^''^ + 6 6';^^'^ + e^ ^^ufJR ^ ^(^3)^ 






(5.1b) 



(5.1c) 



with Cj and h- representing general coefficients in the expansion over e. 



As a first step, we perform the coupling constant renor- 
malization in the MS renormalization scheme. Note that 



in the functions V\jy{t) and r^^(e), defined by (4.7), the 
singularities are contained in 



while the remaining artifacts of dimensional regularizati- 
on can be found in 



r(i-^) 

r(l-2e) 



(47r)'^ = 1 + e(-7 + ln(47r)) + 0{/) (5.2b) 



r(6)r(l-6) = 



1 



sin TT e e 6 

-r(-6)r(i + .) 



■e + 0(e3) 



(5.2a) 



(and similarly for \^\jy and F^^ functions (4Ji3)). By 
expanding F functions over e, in relation (5.1), an un- 



necessary complication of keeping track of various 7, tt^, 
and ln47r terms, would be introduced. Instead, we make 
use of the freedom in defining the MS scheme beyond 
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©(e"), which is explained in detail in Appendix y (along 
with some other conventions and "misconventions" ) , and 
define the bare coupling constant as in terms of the run- 
ning couphng asiy^'^) by 



as 



f^R 



eV'j^lie)] as{f^'n){l 



47r 



1 



-/3o- 



(5.3) 



where /3o — 11 — 2/3 n/. The inclusion of the factor 

tT\jy{e) in ( |5.3| ) turns out to be very suitable for this 
type of calculation, in which both UV and IR singu- 
larities are regularized by the dimensional regularization 
method. The elegance and ad vantage introduced in the 
calculation by the choice (5.3) becomes clear when one 



notes that 



(e) 



n(l) 
^ UV 

rS(e) 



.(0) 



r^'^(6) + o(6^) 



(5.4a) 
(5.4b) 



So, one can see that the presence of the factor er^y(e) in 
( ^.3[ ) is natural in the sense that it contains the combina- 
tions of r's that naturally emerge in this calculation, and 
leads to their cancellation without expanding the whole 
result over e. That is in contrast to "artificial" choices 
like exp (e(— 7 -I- ln47r)) and (A-kY /T{\ — e) found in the 
literature (for example, pTl E8| and [E9l, respectively). 



J 



By substituting (5.4) into (5.1), and performing the coupling constant renormalization according to (5.3), one 
obtains 



M^M 



where 



(0) ^ ^!£04l_^(i) ^ «s(Mfl) 
47r (47r)2 



^"/ 






(5.5a) 




(5.5b) 



(5.5c) 



Note that the only artifact of dimensional regularization 
we are left with is the dimensional parameter e. The 
result (5.5) is given in a simple and compact form in 
which all terms in the expansion over e are still retained. 
Also, the distinction between the singularities of UV and 
IR origin is still preserved. 



If the coefficients of the 1/e and 1/e^ poles of UV origin 
are different from zero, the additional renormalization 
should be performed, as in the case of the multiplica- 
tively renormalizable composite operator from (j){u, t) 
(3.6). The UV singularities are then factorized in the 
renormalization constant ZM,ren- After all UV diver- 
gences are properly renormalized, the remaining l/(— e) 
and l/(— e^) collinear poles should be, at some factoriza- 
tion scale, factorized in Zm.coI- 



B. Renormalization of the hard-scattering 
amplitude T 



We shall now apply the results of the preceding sub- 
section to the renormalization of our re sult s for t he hard - 
scatte ring amplitude T given by Eqs. (|0|) , (|45|), ( [4.10[ ) 
and, (4.21). By comparing (4.1) and (5.1a), we can iden- 
tify M = T{u,Q^)/{NT/Q^)Tnd Al(^r«(u). 



.(0) 



UV(IB) 



appearing in (5.1b) are 
determined from (4.10|) by expanding the coefhcients of 



^UY iji{e) over e, while the n/-proportional NNLO coefh- 



nf,UV(IR) 



in (5.1c) are obtained from (4.21) by ex- 



cients h 

J 1 r 

panding the coefficients of ^jjv^^y^uv IR. 



(e) over e. The 
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TABLE IV: The NLO coefficients af^'^^^ defined in (^b|) and determined for M^^^ = T'^^\u) (EtH). The n/-proportional 



NNLO coefficients b 



nf,UV{IR) 



defined in (5.1c) and determined for 7\/((2,n^) _ ^(2,™^)^^^ ^^ jl) 






uv 

0.2 



Cf- 



1 — It V 2 u 



+ - hi(l -u)] +{u^ 1-u) 



Cf 



T^Q-G4)i'^(^-")-2l^-'(^--))+^^^^-^^ 



IR 



0-2 



Cf- (3 + 21n(l-u)) + (u^ 1 - u) 

1 — It 



Cf 



19 
1-u V 2 



7) + — — -] hi(l - u) - ln'(l - u) ) + (u ^ 1 - u) 



1 /37 / 405-29 20(5 -15\, ,, , f 85 - 7 85 - 1 



In^(l-M) 



^1 



^„„c/v 



Cf —!— f 1 + - hi(l - u) ) + (u ^ 1 - u) 
1 — It V u 



Cf 



11 
1-u \Y 



2-^1 ln(l -u)-- ln^(l -u) + -U2{u) ] + {u ^ 1 - u) 



25 



u 



,nfJR 



, rif ,IR 



^„„IR 



Cf (3 + 21n(l-u)) + (w^ 1 ~ u) 

1 — It 



Cf 



1 



39 



l-ii V 2 
+ {u -^ 1-u) 

311 



24(5 - 34 24(5 - 18 



+ ■ 



3u 



ln(l - u) - 21n^(l - u) + 2U2{u) 



Cf 

+ 



, 408(5-347 408(5- 270 \ , ,, 

+ ^ In(l-u) 

1 - u V 4 V 9 9u ' ^ ' 



24(5-34 245 -18\ 2 ^ 3 

In (1 - m) + - In (1 - u) 



3 3u 

, 24(5-52 245 -18\,.,, ,, . , , ,„ ,., , 
+ I h 1 Li2(it) - 2Li3(u) + 4S'i,2(u) + (tt ^ 1 - u) 



3 



3u 



special cases of the generalized Nielsen polylogarithms 

ln(l - x) 



Li2(u) = Si,i{u) = 
Li3(u) = 'S'2,i('") = 



dx ■ 
2; 



S'l,2(u) 



1 



da; 



In^(l-x) 



(5.6) 



/o 2; 

appear in these results, an(i the useful identity is 
Si,2{u) = -Li3(l - u) - ln(l - M)Li2(M) 

-i ln(u) ln2(l -u) + '^ ln(l - u) + ^(3) . 

(5.7) 

So, after the coupling constant renormalization 
has been performed, the hard-scattering amplitude 



T / [Nt / Q'^)tdkBS the form given by (5.5), where the LO 
contribution T^^^ is given by (4.5) and the coefficients 



a ■ and 6" ^ ' are listed in Table |IV|. As ex- 

pected, the coefficients of the UV poles in (|5.5| ) vanish, 
since all UV singularities get removed by t he c oupling 
constant renormalization. According to Eq. (3.1), i.e.. 



T{u,Q^) = THix,Q^,^%) (g) ZT,coiix,u;fj,%) , 



the remaining singularities of the coUinear type factorize 
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at the factorization scale /z^ in 
ZT,coiix,u;n%) 

47r \fij. 



~IR 



tto (x,m) 



"|(a*?j) /a^I,' ^' 



(47r)2 V/4 



1 



3 * 






^i^^ -«(-(.,.) 



(5.8) 



where the coefScients al^{x,u) and h^^' {x,u) satisfy 
the relations 

With the help of 



as(Mfl) 
= (H) a5(A^^) 






47r 



(5.9) 



one can easily demonstrate that Zt,coi(x^ u; fi'jp) is indeed 
independent of the hard-scattering renormalization scale 



^^R■ 



ZT^coiix,u;fj.%) 
= S{x — u) 



asifJ-l) 1 ~ 



47r 



a-o ix,u) 



(47r)2 ^e- 



3 ^ 



'TrifJR, 



b7' {X,U) jj^ 



2 



a-o {x,u) 



'-^^-4^-a{-ix,u) 



(5.10) 



After factorizing the coUincar singularities from (5.5) by (5.8), and taking into account that b^^'^"' — a^'^ , we obtain 



,nf,IR _ jji 



TH{x,Q\y.l 



T(0)(x)4 
o^Wr. 



as(M|) 



47r 



(ar(^)-<(^))"ar(-)ln^ + 0(.) 



M> 



Q^ 



+ 



(47r)2 



,n!..UV 



{x) 



.rif.IR 



ix) 



+ (ar(^)-«(^(x)-an.-)ln^)ln^ 



(b^'^^x) - 2a'Ax)) 



In 



Mf 



-'^ IR(^\^„'2' ^-F 



Q2 ■ 2<(^)^^ 



(ar(^)-4^(^)) 



0(6) 



(5.11) 



where the 0(e) terms can now be safely neglected (notice that we have kept all e" terms till the end of calculation). 

Finally, having evaluated all the necessary terms, we summarize our result for the hard-scattering amplitude 
Th{x, Q'^.^^p) in the form 



(0), 



TH{x,Q\^i'p) = T'^>{x,Q 



«s(A^fl) rp(l) 



All 



r)/^(x,Q^/i^) 



(47r) 



-3"/] ^H^"^'(a^>'^^MLAiF) 



(5.12) 



where 



Nr 



Wi^.Q^) - -4tA^°\x)^{x^I-x), 



T^j {x,Q ,fJ-p 



Nt 



(x ^* \ — x) 



-1^ H \^tQ i/^i?,iMF) ~ ~?Y2 



(A«(.)-A«(x)ln^)+(. 

+ [X'^l-x). 



We have introduced the functions A^'^'^ and ^^o;, which are given by 

1 — a; 



and 



1 — a; 



l-x 



9- (8(5-7) ln(l-x)+ln (1-x) 



1 -a; 



457 / (485 -95) (-165 +19) 
"2r^ I 18 ^ 6x 



ln(l - x) 



f-i)'°^"-->4'°'<'-^' 



+ 



while 



yl 



4ol(^) 



(2:"/) 



Cf 



l-x 



(125-26) (45-4) 



(3 + 21n(l-a;)) 



Li2(x) + Li3(a;) - 25i,2(a;) 



r'(^ - c,- 



(245-8) (-85 + 6) 



ln(l - x) + 2Li2(a;) 
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(5.13a) 
(5.13b) 



(5.13c) 



(5.14a) 



(5.14b) 



(5.14c) 

(5.15a) 
(5.15b) 



The colhnearly singular terms removed from (5.5) by 
(^.8[) correspond to 



— |yl(°)(u) ' "^^^^) ^ /1(l)^.^ I "s(m|) 1 






Q 



All 



(47r) 



3-/ 



1 



,(1) 



1 



^^ioK-)+^^<r^(-) 



l(2,n/). 



C. Renormalization of the perturbatively 
calculable DA part 

1. General analysis 

Next, we proceed to renormalize the (/)(m, t) following 
t he p rocedu re ou tlined in subsection V A . By comparing 
(^ with (glU) we identify M = (l)iu,t) and M'^'^ = 
(j)^^^ {u, t), while the Q^ scale corresponds to the scale Q^. 
The coefficients 



+ (ti ^ 1 - li) . 



(5.16) 



-UV 



,nf,UV 



The functions A^^^ ( 5.15 ), whi ch ap pear in ( 5.16 ) and as 
coefficients of In^liip/Q^) in (5.13|), are obviously con- 



,nfJR _ ,nf I -, 

5/ =6/ (u, t) 



(5.17) 



nected to collinear singularities of the hard-scattering 
amplitude T . 



are determined from Eqs. (4.28) and (4.29) by ex- 
panding the coefficients of ^ \jy i^jy^ over e, and from 
Eqs. (4.33) and (4.34) by expanding the coefficients of 
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.(0) 



-(1) 



uvi^)^uviR(^) '^^^^ ^' respectively. Although in this 
work only the n^-proportional part of the NNLO con- 

tribution </)„j: {u, t) has been determined, our symbolical 
analysis can be extended to the whole NNLO contribu- 
tion 0^2^ {u, t). In this case the general coefficients of the 
order a| term 4>^'^\u,t) = A^^^^ 



appear. After the coupling constant renormalization, the 
distribution amplitude 4>{u, i) = A^ is given by the ex- 
pression (5.5) with the renormalization scale denoted by 
jlj^. The complete order a|(/i|j) coefficient jd^^^ reads 



I.UV 



lIR 



bi{u,t) 



(5.18) 



AJ(2) ^ 




4|-*(§l «. 



Al 



"¥-'h|' - 



^-*(g|) «2]+0'l.') 



-2 \ 2e 



Q' 



(5.19) 



r 



where Eqs. ( 5.17 - ^^18 ) have already been taken into ac- 
count. 

As denoted in (U|) 

4){u,t) 

the remaining UV singularities are multiplicatively renor- 
malizable and factorize in the renormalization constant 

^0,re«('", v; jl%) given by 



1 



47r 



■oo 



Qs(Mfl) 1 



(47r)2 £2 



6o 



/?oao 



bi 



/3o ai - flo ai 

(5.20) 



with 



bo-Poao-al^O (5.21) 

(i.e., &o(a;,2/)-/3oao(a;,2;)-ao(x,u)(8)ao(M, 2/) = 0) emerg- 
ing as the condition of multiplicative renormalizability. 
As for the coUinear singularities, they factorize at the 
factorization scale ^^ in Z^^coi{s, t; /Iq) given by 



= 1 

X 



as(Mo) 1 „ , a|(A^o) 1 



47r 



(4^) 



— ~ Poao - alj + e f — - /3o ai - ao ai 



(5.22) 



and (5.2C-5.22), the function 



Finally, based on Eqs. (3. 

(j)viv, s; /i|j, Mo) is obtained. It is free of singularities, and 

after the e ^ limit is taken, it takes the form 



= 1 



47r /i;^ 


1 «l(Al?,) 
(47r)2 






^0 ,2 /ifl 
2 /i2 


+ (foi- 


- 2/3o ai - 2ao 


aij 


In'^f 



+ ■ 



(5.23) 



Note that the auxiliary scale Q^ has disappeared after 
renormalization and factorization of coUinear singulari- 
ties. We can make a distinction between the scale jlj^ ^ 
introduced by the coupling constant renormalization and 
the scale /i|j 2 ^t which the remaining UV singularities 
are factorized in the renormalization constant Z^^ren- It 
can be easily shown that the scale /x^ i vanishes from 
the end results, i.e., that Z^^ren and (j)v depend only on 
the scale /i^ 2 ■ Hence, fcR = jl\ 2 and /Xg are the only 
relevant scales. Also, note that Z^^coi{t^'^) — Z7\^^{^'^), 
which is expected, since, in dimensional regularization, 
(^ = 1 when the distinction between UV and coUinear 
singularities is abandoned. 



2. Remarks on the evolutional part of the DA 

As explained in Sec. II, the function (j){u, t) represents 
a perturbatively calculable part of the unrcnormalized 
pion distribution amplitude $(u). By taking into ac- 
count Eqs. (|3.7D and (|3.^) the distribution $(u) can be 
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expressed by ( |3.9| ) 

where $(s, /ig) represents the pion distribution anipHtude 
determined at the scale /ig. Its evolution to the scale /i|, 
is determined by 0y(D, s;/2|j, /ig) and given by (3.11): 

The evolut ion potential V defined in (p]^) can be ob- 
tained from (3.5) 



V 



-z 



-1 

tp,ren 



d 



t^RQ-rZ. 



(f>-,ren 



using (5.20), and it reads 






(47r)2 



v. 



(5.24a) 



where 



Vi = an 

V2 = hi — 2/3g ai — 2ao ai , 



(5.24b) 



By noting that 



47r ^xl 



/3o 






0{a%) (5.25) 



and by e mploy ing the multiplicative r enorma lizability 
condition (5.21) as well a s the results (5.24b), the as 
expansion of (/)y given in ( 5.23| ) can be reorganized and 
written in the form 






(5.26a) 



where 



6^° 



and 



2^1 1 



47r 



Vi 



«i(MMi^.^^(^z 



(47r)2 



iNLO 
9y 



M^ 2 

'a5(MlfJ 
47r 



/Jo'^i) 



In4v2 
Mo 



(5.26b) 



(5.26c) 



denote the L O and NLO part, respectively. By substi- 
tuting ( ^.26a| ) into ( 3.11 ) one obtains 



(^Sifiji) ,NLO 






(/.{)!^^(w,s;/i2j,^g)(^$(s,/i2) 



47r 
iO/ , --2 



«5(Al) ^NLO^^~^2^^ 



An 



(5.27) 



As it is seen from (5.26b) and (5.26c), the results of the 
two-loop calculation correspond to the first terms of the 
LO and NLO contributions to the (pv function. 

The complete LO and NLO behavior of 0y(w,s;/i^) 
and, consequently, of $(u, /2 p) c an be determined by so lv- 
ing the evolution equation ( |2.5| ), or equivalently ( p.l2| ) 

~2 c* , ~2 2\ 

= V^(w,s',M^) «> 0v(s',s,/ifl,Mo)• 
The LO result is of the form 



L.LO 



V (">s;Ai?J = 



n=0 



iV„, 



'-Cl''{2v-l) 



y.Cl'\2s-l) 



as(Aio) 



7/9o 



(5.28) 



where A^„ = (n+ l)(n + 2)/(4(2n + 3)), and Cf/^(2a;- 1) 
are the Gegenbaucr polynomials representing the eigen- 
functions of the LO kernel Vi with the corresponding 
eigenvalues 



7i 



(0) 



Cf 



n+l 



(n+l)(n + 2) 



^H) 



(5.29) 



One can show the agreement between the comple te LO 
prediction given above and the expansion (5.26b). The 
complete formal solution of the NLO evolution equation 
was obtained in p3] making use of conformal constraints. 



and the form of ( 
listed in |1. 



r'-' can be extracted from the results 



3. Analytical results up to n/ -proportional NNLO terms 
(obtained using the naive-ys scheme) 

After this lengthy general analysis we now turn to 
displaying the results. Using the multiplicative renor- 
malizability condition (5.21) and the notation (5.24b), 
the renormalization constant Z0^ren(M, w; /i^) from Eq. 
( |3.20D is expressed by 



5{u — v) + 



"s(Afl,) 1 . . ^ , a|(Afl) 1 



47r 



■Vi{u,v) + 



(47r)2 



2 \ -Viju.v) , V'''{u,v) 
-:;nf ] { +e 



+ ■ 



(5.30) 



Here we list only the relevant combinations of a, , &j co- 
efficients: 



Vi{u,t) = ao{u,t) 



^ 2Cp I- 



t — u 



e{t-u) + 
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V^f {u,t) ^ b^;^ {u,t) - 2ai{u,t) 

1 \ A5 
3 




In - 1 e{t- u) 

(5.31b) 



Our results confirm the well-known form of the one-loop 
kernel Vi \ J and the two-loop nj^-proportional kernel ^2"^ 
p9| , pO| pi I . For later use we als o specify the convolution 
of the functions given in ( ^.31 ) with the frequently en- 
countered 1/(1 — x) term: 



(g) Vi{x,u) = Cf 



l-u 



(3 + 21n(l-u)) 



and 



X 

Cf 



V,"'{x,u) 



16 

y 



ln(l- 



2Li2(M) . 

(5.32b) 



4- The results obtained using the HV scheme 

In the following we present the results obtained using 
the HV scheme. 

The function (j) calculated in the naive-75 scheme and 
the function (f)^^ obtained in the HV scheme are related 
by 






(5.33) 



The factors Z^jy jjy and Z^y ^^j remove the "spurious" 
anomalous terms introduced by the presence of dimen- 
sionally regulated UV and coUinear singularities, respec- 
tively. 

The renormalization of the (j)^^ function proceeds 
analogously to the renormalization of (f> described in pre- 
ceding subsection. While 



-HV 



lHV 



the coefficients Oj, bi for i > 1 get replaced by 



6f^ 



A6,, 



(5.34) 



(5.35) 



and, according to ( |3.8| ), the UV and coUinear singularities 
are factorized 



IHV _ yHV iHV yHV 



(5.36) 



By c omparing Z^^„, Z^^,, and 0^^ with the re- 
sults (5.20-5.23) obtained using the naive-75 scheme, Eq. 
(^.36|) takes the form 



IHV 



{■^4>,ren '^HV,UV 



)(^ 



rfin 
■HV.,UV 



(t>V Z 



fin 
HV,col 



\ZhV,co1^4>,coi} 



(5.37) 



where 



rjdiv czdiv 

^HV = ^HV,UV 



alifi 



s\^^Rl 



(Att) 



A61 



/SoAai — ao Aai 



(5.38) 



(5.32a) and 



Z 



fin 
HV 



czfin 
^HV,UV 



47r 
A61 



(Aai + 0{e 



«l(Afl) 



(47r)2 
/3oAa2 - flo Aa2 - ai Aai + 0(e) 



(5.39) 



while Z'^y'll'lifj.'^) = (Z'^v.uv) ^(m^)- The condition of 
multiplicative renormalizability of the "spurious" anoma- 
lous terms introduced by the HV scheme reads 



A61 - /3oAai - 2aoAai = . 



(5.40) 



Finally, we list the results obtained by substituting 



'(1) 



42) 



K'hv and Kl,„y 



3WM) in place of K^^^ and i^^^) 



in Eqs. (4.22-03). The combin ations of the Aa j and 
A6,; coefficients, which appear in ( 5.38 - 5!4^ ) after ( 5.40| ) 
is taken into account, read 



Aai{u,t) 



2Cf iA-e{t-u) + 




(5.41a) 



Ab'l'(u,t)-2Aa2{u,t) 



(5.41b) 

The comp lete renormalization constant Zhv = Zhv,uv 
from Eq. ( p.33| ) is then given by 

ZHv{u,v;ilJi) = Zfj'^iu^w;^!^^) (g) zl^y{w,v;il\) . 

(5.42) 

By utilizing the condition of m ultip licative renormaliz- 
ability of the "spurious" terms ( 5.4C ) , it takes the form 
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ZHviu,v;ill) = J(u-t;) + ^4r^(^«iK«) + 0(e) 



(47r)2 



47r 
2 
"3' 



1 -Aai(u,w) A6"^(m,u) -2Aa2(M,u) 



r 



0(.) 



(5.43) 



It is interesting to note that by using the "reduction" 
formulas, which relate the exclusive ( "nonforward" ) and 
inclusive ( "forward" ) kernels |2l[| , the agreement between 
the renormalization constant 2^/ given above and an 
analogous "HV" renormalization constant for the longi- 
tudinal spin structure function gi pl[ is established (up 
to 0(e) terms). 

For later use we s pecify the convolution of functions 
defined in Eq. ( 5.41 ) with the 1/(1 — x) term: 



and 



Aai(x,u) = Cf 



1 



1-u , 
-8 ln(l 



(5.44a) 



1-x 



(Ab"' {x,u) - 2Aa2{x,u)) 



16 1 — u , ,, , „l — u^ . , , 
— ln(l - m) - 8 Li2(u) 



1-M 



(5.44b) 



VI. THE EXPRESSION FOR THE PION 

TRANSITION FORM FACTOR UP TO 

n/-PROPORTIONAL NNLO TERMS 

We now combine the results for the hard-scattering 
amplitude and the DA obtained in the preceding sections. 
After resolving the 75 problem and discussing the depen- 
dence of the prediction on the factorization scale /i|, , we 
finally present the expression for the pion transition form 
factor up to 7i/-proportional NNLO terms. 



A. Resolving the 75 ambiguity 



In subsections |VI3| and |VC| , we have presented the re- 
sults of the perturbative treatment of the hard-scattering 
amplitude T(u, Q^) and the distribution amplitude $(?i) 
respectively. Along the lines outlined in Sec. 



Ill 



proceed to combine these results to obtain the finite and 
75-scheme independent expression for the pion transition 
form factor F^-^.,^{Cp), up to the n/-proportional NNLO 
contributions. 

The lack of the ambiguity in the DA results along with 
the fact that the prediction for the pion transition form 



factor should not depend on the choice of the 75-scheme 
make it possible to resolve the ambiguity of the 75 treat- 
ment in the hard-scattering calculation. 



1. Naive "/5-scheme 

The appearance of two 75 matrices imposes the use of 
the naive-75 scheme in the DA calc ulation and the cor- 
responding results are presented in (5.2(;-5.31). The 75 
matrix present in the hard-scattering amplitude can also 
be treated in the naive-75 scheme in which case a num- 
ber of results emerge. After the Ward identities of QED 
are taken into account, the remai ning ambi guity in the 
hard-scattering amplitude result (5.12- 5.1q ) is parame- 
terized by t he par ameter S (as explained in Appendix H 
and Subsec. |IVA| ) 

Matching Eqs. (5.15) and (5.32) one observes that 



1 



1-x 
1 



l~x 



^Vi{x,u) = A[1]{u) (6.1a) 

,V,"'{x,u) = Af/^^(«)|,=i. (6.1b) 



If these relations are ta ken i nto account in Eq. (5.16) 



then a comparison with (5.30|), for fij^ = fij^ (i.e., for the 



DA $(m) renormalized at the fip scale), gives 



'-JT,col\S=l 



j),ren ' 



(6.2) 



i.e., the relation ( 3.14 ) is satisfied and the singularities 
in (3.13) cancel for the 5 parameter taking the value 1. 
Hence, we obtain 



^T,co. 



i{x,u\y?p) 



-JT,coi 



.i{x,u;n]p)\s=i 



and, consequently. 



Th{x, Q^, ^I) = Th{x, Q^, ti\)\s=i 



(6.3) 



(6.4) 



where Eq. (|6.2| 

Zt.c oi(x,u\h']p)\s 



together with ( ^.30|) , deter mines 
while taking (5 = 1 in Eqs. ( |5.12| - 
5.151) gives Th{x, Q"^, 1^^)16=1- Hereby, we have confirmed 
the 75 prescription employed in pOf (see Appendix H for 
further details). 

In preceding consideration, we have resolved the 75 
ambiguity of the hard-scattering prediction by adopting 
the naive 75-scheme and by using the unambiguous DA 
results (along with the QED Ward identities) to single 
out the correct prediction. 
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2. HV scheme 



while 



Let us now present the calculation performed in the HV 
scheme. Continuing to D dimensions and adopting the 
HV scheme leads to unique results, but the "spurious" 
anomalous terms, which violate chiral symmetry, appear 
and additional renormalization is required, both for the 
DA and the hard-scattering amplitude. 

The results for the hard-scattering amplitude obt ained 
in th e HV scheme correspond to the 6 — choice in (5.12- 
Slel) , and the notation T"^ , T^^ , Z^^^^^ has been in- 
troduced. The fact that the UV singularities appearing 
in the hard-scattering amplitude get completely renor- 
malized by the coupling constant renormalization indi- 
cates that, contrary to the DA case, only the presence 
of coUinear singularities along with the nonanticommut- 
ing nature of 75 matrix introduces "spurious" anomalous 
terms. 

The corresponding renormalization constant for the 
DA denoted by Zhv = 2f}y Zj^y and displayed in (5.38- 
5.43) has been determined by comparing the "correct" 



results obtained in the naive- 75 with the corresponding 
results obtained using the HV scheme. As a result, one 
finds that the unrenormalized DA in the HV scheme, 
^^^{u'), and the unrenormalized DA in the naive-75 
scheme, $(u), are related by 



$ = Z 



HV 



<i>^^. (6.5) 

Similarly, the renormalization constants Z^jy^en a-nd 
Z¥Yf,n determined in the naive-75 and the HV schemes, 
respectively, are related by 






1 r^HV 

ip.ren ' 



(6.6) 



while the additional finite "HV" renormalization of the 
renormalized distribution ^^^{v',jsj^) calculated in the 
HV scheme, is needed to obtain the renormalized DA 
<i>(ti,/2|j.) free of "spurious" anomalies: 



$ 



2;/™-! (f, 



HV 



(6.7) 



The prediction for the pion transition form factor can- 
not depend on the choice of the scheme and chiral sym- 
metry is restored for the complete result, i.e., 



F^.^^ = T'^^ =TZ„l Zhv ^^ 



(6.8) 



or equivalently 



F . 



Th'^* 



Th [Zl^ 



2-fin 
■-HV 



$* 



Tg^ <P"^* 



(6.9) 



On the basis of Eqs. (|5.44|) and (|5.38| - ^.43| ), along with 
Eqs. (5.12-5Tl6|), it can easily be shown that 



T = T"'' Zhv, 



(6.10) 



rp rpHV -J fin 

-L H — -t w ^ 



^H 



HV ' 



and 



r/HV 'j'div 
^T,col ^HV ' 



(6.11) 



(6.12) 



with Th and Zt,coI being given by Eqs. (6.4) and (|6.3|), 
respectively, and T = Th Zt,coI ■ 

Hence, we have resolved the 75 ambiguity appearing 
in the hard-scattering amplitude calculation by consis- 
tently treating, in either the naive-75 or the HV scheme, 
both the hard-scattering amplitude and the distribution 
amplitude (which, actually, is free of the 75 ambiguity). 



B. Discussing the factorization scale dependence 

After the 75 ambiguity is resolved, and the coUinear 
singularities present in the T{u, Q^) and <&(«) amplitudes 
cancel, we are left with the finite prediction for the pion 
transition form factor 

F^.^^{Q^)^TH{x,Q^,til) ® <^*{x,fil). (6.13) 

The hard-scattering amplitude Th{x, Q^,fJ.^), evaluated 
up to rif- proportional NNLO terms, is given by Eqs. 
( ^■12[j5.15| ) with 6 = 1. 

The distribution amplitude <i>*(x, /i|,) is determined by 
evoluting ^*{x, /Iq) (obtained at the scale /Iq using some 
nonpe rturbative method) to the scale ^^ according to 

$(W, ill) = 4>v{v, S- fll, ^2) ^ $(g^ ^2^ ^ 



i.e 



*(",Mf) 






47r 



i^^O, 



In subsection V C we have analysed in detail the evo- 
lutional part (j)y, and as noted there, the two- loop DA 
calculation explicitly gives only the first few terms of the 
LO and NLO contributions, while the complete LO and 
NLO behavior of t he D A is determined by s olving the 
evolution equation ( p.5\ ), or equivalently (3.12) 

d 



fJ-l' 



= V"(w,s',/i^) 



The dependence of Th{x, Q^ , fip) on the factorization 
scale ^^ can be determined analogously to the / ip d epen- 
dence of the DA. Thus, by differentiating Eq. ( |3.lD with 
respe ct to ^^ and taking into account Eqs. ( |3.5p and 
( |3.14|) one finds that Th{x, Q^,fJ-j?) satisfies the equation 

l^'F-;rT^H{x,Q'^,y^p) 



df^i 



-TH{y,Q^,fj.F) (E> V{v,x-iil), 

(6.14) 
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which, as it is seen, is analogous to the DA evolution 
equation (p!q). Therefore, j ust as in the case of the DA, 
any finite order solution of (6.14) contains the complete 
/i|n dependence of T h{x, Q^, fij^) which is not the case 
with the expansion ( ^.12 ) tr uncat ed at the same order 
(and not taking into account ( 5.25| )). 

The hard-scattering amplitude Tff (x, (5^,/i|,) can be 
written in the factorized form 



TH[x,Q\til) = TH{y,Q^fil = g' 



•<t)v{y,x,Q'^,fj,'jp), 
(6.15) 



with <pv {y , X , Q'^ , ^jl^) containing all its /i|, dependence. 
This can easily be demonstrated if use is made of T^ and 
(j)v de term in ed up to n f-pro portional NNLO terms (see 
Eqs. (^.26|) , ( |5.32[ ), and (|C3b|)). On the other hand, usin, 



Eq. ( 3.12 ) one can show by partial integration that ( 6.15 ) 
i ndee d represents the solution of the evolution equation 

(HI- 

When calculating to any finite order in as, it is inap- 
propriate to convolute the <&(a;, fi %) di stribution obtained 



by so lving the evolution equation ( p.5[ ) (i.e., given by (2.8 
2.13 )) with Th {x, Q ^, /ip) obtained by the truncation of 



the expansion (5.12). Namely, in the latter function, only 
the partial dependence on fi^ is included, in contrast to 
the former. Notice that when the complete dependence 
of Th on /i|n is taken into account even the LO term in 



(3.15) is is^ dependent, in contrast to T^^^ given in (2.4), 
and this leads to ^^ independent LO prediction for the 

pion transition form factor. 

Substituting Eqs. ( |6.15D and ( |3.1l| ) into Eq. ( |6.13| ) 
and taking into account that 



(j)viy, X, Q^,tiF) (g) 0y (x, s, flp, ^l) = 0y (y, s, Q'' 



(6.16) 



one obtains 



F^•7-(Q')-^H(y,Q^Q' 



(6.17) 



The relation (6.16) is valid at every order of perturba- 
tive calculati on (t o NLO this can easily be checked by 
substituting (5.28) and the NLO results of Ref. |23] into 
( |6.16|) ). It represents the resummation of the ln(Q^//iQ) 
logarithms over the intermediate /i|n scale, performed in 
such a way that, first, the logarithms ln(^|,/^Q) originat- 
ing from the perturbative part of the DA are resummed, 
and then the summation of ln((5^//i|,) logarithms from 
the hard-scattering part is performed. Therefore, the 
summations of the /i|n logarithms can be accomplished 
with any choice of jj^p, because the effect in the final pre- 
diction, at every order, is the same as if the complete 
renormalization-group resummation of \n[Q'^ / fi^) loga- 
rithms has been performed. 



Consequently, the Fj*^,^ prediction (2.2) (as well as 
the prediction for any other exclusive quantity obtained 
in the standard hard-scattering picture) is independent 
of the factorization scale ^p at every order in as, pro- 
vided both Th and <I> are consistently treated regarding 
the /i|n dependence. The intermediate scale at which the 
short- and long-distance dynamics separate, the factor- 
ization scale /x|i, disappears from the final prediction at 
every order in as and therefore does not introduce any 
theoretical uncertainty into the PQCD calculation for ex- 
clusive processes. 

We would like to point out here that by adopting 
the common choice /i|n — Q^, one avoids the need for 
the resummation of the ln((5^//i|') logarithms in the 
hard-scattering part, making the calculation simpler and 
hence, for practical purposes, the preferable form of 
Fj'^yriQ^) is given by 



F^,^^iQ^)^TH{x,Q^,Q^) ® $*(x,g2) 



(6.18) 



C. Presenting the final results 



Finally, we summarize. Taking into account Eqs. ( 5.12 - 5.15 ) and the results of the preceding subsections, the hard- 



scattering amplitude Th{x,Q ,fip = Q ), free of all effects of coUinear singularities (the terms containing functions 
^coi ^'^" Mf/Q^ factorized in (t>v), takes the form 



THix,Q^,fi%^Q^) = tI^\x,Q 

(47r)2 



"^(^^)TW(x,Q^M^ = g^) 



4tt 



-3"/ 



.(2 



^r'^^(x,Q^M^ = Q^M^) 



(6.19) 



where 



TPix,Q') = 



tI,'\x,Q',^i%^Q') 



T\^^'''\^.Q\^4 = Q^^^l) 



Q 



2 A(°\x) + {x^l-x), 



Nt 
12" 



Q 






mI 



+ (x ^ 1 - .t) 
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(6.20a) 
(6.20b) 
(6.20c) 



and 



A^^'^x) = 


l-x' 


A^^\x) -- 


1 ~ X 


A(2,«/)(a;) , 


I — X 



-9 - ln(l - x) + 111^(1 - x) 



457 /47 1 \, ,, ^ A3 1 \, 2,, 
iln^(l - x) - yLi2(x) + Li3(a;) - 25i,2(a;) 



(6.21a) 
(6.21b) 



(6.21c) 



To obtain the distribution amphtude $(a;, Q^) one evo- 
lutes $(a;,/iQ), determined using some nonperturbative 
methods at the scale /Iq, to the scale Q^ according to 
(p[]2l3|). 

By substituting Eqs . ( |6.19| ), (U), and ( |2.10| ) into 
( |6.18 ) and taking ( C3b ) into account the pion transition 
form factor Fy*jT^ expressed as a perturbative series in 
as(/x|,) reads 






„2/',,2 



"s(AiflJ 5,(1) 



(47r)2 



47r 
2 



^7*77r W 



,^^-r(Q^M^) + ■ 



(6.22) 



where 



-^7»77r W 
-^7*77r W 



_ ^(0) 



) = T),"'(x,Q")®$^^(:E,g'), (6.23a) 



.(0 



^Q^)(»$^^(a;,Q^) 



+rJ^'^(a:,Q^)®$""^(a;,Q^), 



(6.23b) 



p(2,nj,).„2 ,2 N, 



Ti,''"^^(2:,Q^Q^A^|^)®$^«(2:,Q2) 
+ri,°Ha;,Q2)^$A^LO(^g2)ij^^^ 

(6.23c) 



As it is seen, we are left with one expansion parame- 
ter: as(/^|f). The scale /i|j actually represents the renor- 
malization scale of the complete perturbatively calculable 



part of the pion transition form factor (3.17), i.e., of 



THiy,Q^tit,) = THiy,Q,Q')®^v{y,s,Q',f,'o). 

(6.24) 

Although, the physical pion transition form factor 
Fj»jTr{Q^) does not depend on the choice of the renor- 
malization scale /k^j, when calculating to any finite order 
a residual dependence on the /^^ scale remains. 



VII. NUMERICAL PREDICTIONS 

A. Fixing the renormalization scale according to 
the BLM procedure 

The dependence of finite order predictions on the 
renormalization scale introduces a theoretical uncer- 
tainty in their interpretation (see ||] for a detailed dis- 
cussion), which is especially evident in calculations to 
lowest order in order as- It would be advantageous to 
optimize the scale choice according to some sensible cri- 
teria. The BLM procedure 11^, O, |l5| offers such cri- 
teria. The essence of the BLM procedure is that all 
vacuum-polarization effects (giuon vacuum polarization 
contributions, analogous to QED, as well as, quark vac- 
uum polarization and vertex corrections) from the QCD 
/3 function are incorporated into the running coupling 
constant rather than into the coefficients of the pertur- 
bative expansion. In practice, this amounts to computing 
quark-loop insertions in the diagrams of that order (since 
/3o — 11 — 2/3nf,- ■ ■ ) and setting the scale by demanding 
that the coefficients of the perturbative expansion are nf 
independent. 
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Hence, according to the BLM scale setting prescrip- 
tion, the renorniaUzation scale for the pion transition 
form factor entering at the NLO is determined from the 
NNLO n^-proportional terms, and is fixed by the require- 
ment 






Note that in the present calculation the effective nature of 
I 



/i|j. has been implicitly assumed, i.e., /i|j, has been treated 
as independent of the momentum fractions x through- 
out the paper; otherwise the factorization of singulari- 
ties would take a cumbersome form (if menagable at all). 
Apart from that we would be faced with the problem of 
a clear separation of the short- and long- distance effects 

i5i. 



Therefore, the only consistent way to asses the BLM scale is to solve Eq. (7.1), resulting in some mean value BLM 
scale 



where 



asLMiQ^) = exp 



M_R = t^BLM = o,blm{Q ) Q , 

A(^."/)(z)(^$^o(x,Q") 

yl(l)(.T) ® $L0(X,Q2) + A(")(x) ® $ArLo(x,Q2) 



(7.2a) 



(7.2b) 



As it is seen, the scale /i^^^/ depends on the specific form of the distribution amplitude. 

In Section II the nonperturbative input, i.e., the distribution amplitude (jjix,^"^) — ^( x, p.a) /Nd, determined at the 
scale /Zq, has been presented in the fo rm of an expansion over Gegenbauer polynomials ( 2.1C ). The evolution to the 
scale fij^ has been desc ribed by Eqs. (2.11-2.13). Retaining only the first three terms in the general expansion of the 
pion DA given in Eq. ( 2.10 ) 



3/2/ 



-i3/2. 



(x,Aio) =6x{l-x) l + B^C^' {2x-l) + BiCl' {2x-l 



the LO and NLO contributions to Q^F^^^T^iQ^) take the form 



Q'fP^JQ') = 2C^U[3{l + B^^iQ^) + BiOiQ'))] 






2(~^ttJtt 



-20 + ^s,^°(g^) 



10487 , 
225 



Bi^iQ^^+iY^'B^'^'^Q" 



k=2 



(7.3) 



(7.4) 



(7.5) 



while the Uf proportional NNLO contribution amounts to 



fF^l-;^\Q^ 



^C^U [(-43.47 + 78.47S^°(g2) -|- 197.165Bf '^(Q^)) 



-20+^-^B^^iQ') 



225 



-srm)in^ 



3 5: 'i?r°(Q^)j In ^ 



(7.6) 



Therefo re, o n the basis of Eqs. (7.1) and (7.6) one finds that the BLM scale for the pion transition form factor is 
given by ( [7.2a ) with 



/ 



asLMiQ'^) = exp 



\ 



-43.47 -t- 78.47S|'0(Q2) + 197.165BfO(Q2) 



-20 



295 

Is" 



B^'^m + ^Sr (Q^) + 3^ 'Sf ^°(Q^) 



(7.7) 



fc=2 



Expressions (6.22) and (7.4-7.6), representing the com- plctc NLO prediction for the pion transition form factor. 
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together with the expressions (7.2a) and (7.7), specify- 
ing the corresponding BLM scale, are vahd for arbitrary 
distribution amplitude (with the evolutional effects in- 
cluded) and represent the main results of this paper. 



B. Numerical predictions in the MS and av 
schemes 

Based on the general expressions derived in a preceding 
subsection, we now proceed to obtain numerical predic- 
tions for the pion transition form factor using two specific 
distribution amplitudes: the asymptotic DA and the CZ 
distribution amplitude. 

There is increasing theoretical evidence coming from 
different calculations ||3^, ^, |3j, Q that the low energy 
pion distribution amplitude docs not differ much from its 
asymptotic form(/)Qs (a;) = 6a; (1 — a;) ( whi ch represents the 



solution of the evolution equation ( |2.5| ) for /zj, -^ oo). 
The distribution 0(a;,/io) = 4>as{x) is characterized by 
the fact that at the LO it has no evolution, while the 
NLO evolutional effects are tiny [B|, and for the purpose 
of this calculation these effects can safely be neglected. 

The expression for the pion tr ansiti on for m facto r 
Q'^F^.^^iQ^), based on the Eqs. (|^ and (0-^, 
and corresponding to the asymptotic distribution then 
reads 



^ 7*7"^ V^ / 



(47r)2 



47r 



(-20) 



-3-/ 



^43.47 -20 In ^ 



{7- 



The n/-proportional NNLO contribution determines the 
value of the BLM scale 



f^R ^ [t^BLM) 



0.114Q2 « ^ 
9 



(7.9) 



One notes that this scale is considerably softer than the 
total momentum transfer Q^, which is consistent with 
partitioning of Q^ among the pion constituents. It should 
be pointed out, however, that in the MS scheme the 
BLM scal e do es not reflect the mean gluon momenta. 
Based on (7.8), the NLO prediction amounts to 



Q2^^.^^(g2) ^ 0.185 



"s(/4.) 



(-0.309) 



(7.10) 



This prediction obtained with the asymptotic distribu- 
tion and the MS definition of the str ong coupling renor- 
malized at the BLM scale given by (7.9) is displayed in 
Fig. g, along with the CLEO experimental data. For 
comparison also included in Fig. || is the NLO prediction 
obtained by employing the widely used choice /i^ = Q^. 
The usual one- loop formula for the QCD running cou- 
pling constant (C5) has been used with A = A-p-j- = 0.2 
GeV^. 



As it is seen from Fig. |^, the NLO results for the 
pion transition form factor display the following features. 
First, inclusion of the NLO contributions decreases the 
LO prediction. Second, predictions based on the asymp- 
totic distribution are in reasonably good agreement with 
currently available experimental data. 

In comparison with the choice fij^ = Q^, the BLM 
scale choice increases the absolute value of the ratio of 
the NLO to LO prediction by « 11-6% for the val- 
ues of Q^ between 6 and 20 GeV^. As an extension of 
the BLM scale-fixing prescription to all orders in per- 
turbation theory, in [^ all the (/Soa^)" contributions to 
the pion transition form factor were resummed under the 
assumption of " Naive Non-Abelianization" (NNA) . Our 
results cannot be directly compared with [ p9[ , but there 
are indications that the usual BLM scale-setting overes- 
timates the size of higher-order contributions associated 
with the one- loop running coupling pO[ . 

The complete NLO prediction for the pion transition 
form factor corresponding to the end-point concentrated 
CZ distribution amplitude |^, given by Eqs. (7.4-7.t) 
with the coefficients B2 = 2/3 and B4 — 0, is shown in 
Fig. [1^. Owing to the fact that the LO and NLO evolu- 
tional corrections to the CZ distribution arc considerable 
H, they have been taken into account. The BLM scale 
for the CZ distribution is higher than the BLM scale 
for the asymptotic DA, and it varies from (3^/1.84 to 
gV2.37 for Q2 between 2 and 20 GeV^. 

As it is evident from Fig. nO^, the complete NLO pre- 
diction for the pion transition form factor derived from 
the CZ distribution exceeds the experimental data signif- 
icantly. This result can be considered as a serious failure 
of the CZ distribution amplitude. Comparing Figs. O 
and 10, one observes that the difference between our re- 



sults for the asymptotic and CZ distribution amplitudes 
is sufficiently large for an unambiguous experimental dis- 
crimination between the two possibilities. Therefore, one 
expects that the pion distribution amplitude is closer to 
the asymptotic form than to the strongly end-point con- 
centrated DA's like CZ 



34 1 



The size of higher-order QCD corrections represents 
the missing ingredient in assessing the validity of the 
perturbative prediction and the convergence of the ex- 
pansion. One hopes that the BLM prescription offers a 
systematic way to choose the renormalization scale and 
minimize higher-order contributions. In order to check 
this for the case of the pion transition form factor, one 
would have to evaluate the complete NNLO contribution, 
which is a very demanding task. Another sensible indica- 
tor of the applicability of the perturbative calculation is 
the size of the expansion parameter as{fJ,Ji)- The rather 
low BLM scale given in (7.9), and consequently the large 



'^sif^'^LAi)' questions the applicability of the perturba- 
tive prediction at experimentally accessible momentum 
transfers. Namely, the NLO predictions obtained in this 
pape r assuming the asymptotic DA and the BLM scale 
(7.9) satisfy the requirement as{lJ,%) < 0.5 for Q^ > 6 
GeV^. It should be pointed out, however, that there is 
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FIG. 9: The LO and NLO predictions for the pion transition 
form factor (scaled with Q ) Q Fj»jn(Q ), obtained using 
the asymptotic DA. The NLO predictions are obtained using 
the BLM scale /x|j = fJ-%LM ~ Q'^/9 (7^) and the commonly 
used choice /i^ — Q^. The experimental data are taken from 



FIG. 11: The LO and NLO predictions for the pion tran- 
sition form factor (scaled with Q ) Q F^'-jniQ ), obtained 
corresponding to th e asymptotic DA in the av scheme with 



fi. 



Q /1. 7 (7.14). The experimental data are taken from 
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FIG. 10: The LO and NLO predictions for the pion tran- 
sition form factor (scaled with Q^) Q^F-ytj^{Q^), obtained 
using the CZ distribution amplitude. The NLO predictions 
are obtained using the BLM scale determined from (ITJ) , and 
the commonly used choice fi^ — Q^. The experimental data 
are taken from fB6|l. 



an intrinsic disa dva ntage in using the AIS running cou- 
pling (given by (|C5|)) as an expansion parameter, since it 
has a simple pole at /i|j — A^. This does not reflect the 
nonpcrturbative behavior of Q.s{p%) for small /i|j, and 
a number of proposals have been suggested for the form 
of the coupling constant in this regime UU, M, but its 
implementation demands caution ||^. For a recent ap- 
plication of 1 42 to the calculation of the pion transition 
form factor, see [[i3[ . 

So far nothing has been said concerning the renor- 
malization scheme dependence of the predictions. It is 
known that the renormalization scheme dependence can 
be avoided by considering relations between physical ob- 
servables, which must be independent of the choice of 



the scheme and scale to any fixed order of perturbation 
theory. This requirement can be expressed in the form 
of "commensurate scale relations" (CSR), in which the 
BLM scale-setting method is used to fix the renormaliza- 
tion scale [|lj|. In practice, a CSR for two observables 
is obtained by relating their respective perturbative pre- 
dictions calculated in, for example, the AIS scheme, and 
then by algebraically eliminating ctjj-g. The choice of the 
BLM scale ensures that the resulting CSR is indepen- 
dent of the choice of the intermediate renormalization 
scheme. Following this approach, in ||l^ the exclusive 
hadronic amplitudes were calculated in the ay scheme, 
in which the effective coupling ay(/j,^) is defined from the 
heavy-quark potential V{fi^). The ay scheme is a nat- 
ural, physically based scheme, which by definition auto- 
matically incorporates vacuum polarization effects into 
the coupling. The fXy scale which then appears in the 
ay coupling reflects the mean virtuality of the exchanged 
gluons. 

If use is made of the scale-flxed relation between the 
couplings a-jjg and ay p^ 



Q^Afs(^ 



-5/Vv) = «v(M^y)(l + 



avip-v) 8Ca 



An 



(7.11) 



then, to the order we are calculating, the prediction for 
the pion transition form factor takes the form 



avif^v) p(i) 

'f'j-iT\^ ) — -f'^'-yTrV^ ) ^ ^ 7*7' 



Kr,AQ') = F\%AQ') 



F\iUQ') + 



with the "V" scale being given by 



,,2 _ p5/3,,2 

^y — e fJ-BLM ■ 



(7.12) 



(7.13) 



Now, on the basis of ( |7.9| ) and (7.13) one finds that the V 
scale, corresponding to the asymptotic distribution and 
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reflecting the mean (NLO) gluon momentum, is 



{^?y^^Qm2hQ'^^ 



(7.14) 



Just for comparison, it is worth mentioning here that the 
corresponding ^y scale for the pion electromagnetic form 
factor amounts to Q'^ /2Q ||, |l|l. 

Furthermore, since ay is an effective running coupling 
defined from the physical observable it must be finite at 
low momenta, and the appropriate parameterization of 
the low-energy region should in principle be included. 
Nevertheless, in the energy region we are interested in, 
the usual one-loop conventional solution of the renormal- 
ization group equation for the QCD coupling (C5) can be 
employed. The numerical NLO prediction for the pion 
transition form factor obtained from Eqs. (7.12-7.14) and 
calculated with Ay = 0.16 GeV^ is depicted in Fig. |l|. 
As can be seen, it is in good agreement with experimental 
data. The LO QCD correction, i.e., the NLO contribu- 
tion, lowers the LO prediction for w 16% for Q^ sa 6 
GeV2, i.e., for ay(Mv) ~ 0.3. 



VIII. SUMMARY AND CONCLUSIONS 

In this paper we have determined the NLO Brodsky- 
Lepage-Mackenzie (BLM) scale and obtained the com- 
plete NLO prediction for the pion transition form factor. 

To determine the NLO BLM scale, a consistent and de- 
tailed calculation of both the hard-scattering and the per- 
turbatively calculable part of the pion distribution am- 
plitude has been performed up to ny-proportional NNLO 
contributions. The calculation has been carried out in the 
Feynman gauge. To control both the UV and coUinear 
divergences the dimensional regularization metho d ha s 
been employed. By combining, according to Eq. (3.13), 
and matching the results for the hard-scattering ampli- 
tude with the corresponding results obtained for the dis- 
tribution amplitude, a proper cancellation of coUinear 
singularities has been established and the 75 ambiguity 
problem (related to the use of dimensional regulariza- 
tion) has been resolved using the naive-75 as well as the 
't Hooft-Veltman (HV) schemes. As a result, the com- 
plete leading-twist NLO prediction for the pion trans ition 
f orm fac tor has been obtained as given by Eqs. ( 6.22 ) and 
( [tIHt^ , and Eqs. ( [7.2a| ) and (^), specifying the cor- 
responding BLM scale. Derived in the MS scheme (for 
which suitable compact form ( |5.3| ) has been adopted), 
these expressions are valid for arbitrary form of the distri- 
bution amplitude (with the evolutional effects taken into 
account), and represent the main result of the present 
paper. 

It has been demonstrated that the prediction for the 
leading-twist perturbative QCD prediction for the pion 
transition form factor is independent of the factorization 
scale ^^p at every order in the strong coupling constant 
as provided both the hard-scattering amplitude and the 



pion distribution amplitude are treated consistently re- 
garding their p.'^p dependence. The factorization scale 
disappears from the final prediction at every order in as 
without introducing any theoretical uncertainty. Con- 
sequently, for practical purposes the simplest and com- 
monly used choice y?p ~ Q^ is justified at the intermedi- 
ate steps of the calculation. 



Based on the general expressions ( 6.22 ) and (7.4-7.6), 
the NLO predictions for the pion transition form factor 
have been obtained using the asymptotic and the CZ dis- 
tribution amplitudes, with the renormalization scales be- 
ing given by the respective BLM scales determined from 
Eqs. ([7.2aD and ( [7.7D . These predictions are shown in 
Figs. ^ and |l^, respectively. By comparing these fig- 
ures, one observes that, while, on one hand, the pre- 
diction derived from the asymptotic distribution is in 
good agreement with the presently available experimen- 
tal data, on the other hand, the prediction obtained as- 
suming the CZ distribution exceeds the data significantly, 
clearly demonstrating the inadequacy of the CZ distribu- 
tion. This is in accordance with the conclusions reached 
in |3|, ^ ^, |3^, according to which the distribution 
amplitude is closer to the asymptotic form than to the 
end-point concentrated distribution of the CZ type. The 
renormalization scale //^ fixed according to the BLM 
scale setting prescription within the MS scheme and 
corresponding to the asymptotic pion distribution am- 
plitude, turns out to be /i|j blm ~ Q^ 1^- Thus, in the 
region of Q^ < 8 GeV^, in which the experimental data 
exist, y^'sLM < 1 GeV^. Consequently, the prediction 
obtained with /i|j 
considered reliable. 



t^BLM cannot, in this region, be 



In addition to the results calculated with the asymp- 
totic distribution and MS renormalization scheme, the 
numerical prediction assuming the same distribution but 
in the ay scheme, with the renormalization scale /ijj = 

My — ^^^^^^'BLM ~ Q^/2j has also been obtained. It is 
displayed in Fig. |Tl| and as seen, is in good agreement 
with experimental data. Due to the fact that the scale 
^y refiects the mean gluon momentum in the NLO di- 
agrams, it is to be expected that the higher-order QCD 
corrections are minimized, so that the leading order QCD 
term gives a good approximation to the complete sum. 
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APPENDIX A: 75 PROBLEM 
1. General remarks 

When using dimensional regularization, one runs into 
trouble with quantities that have the well-defined prop- 
erties only in D — A space-time dimensions, that is, with 
the Levi-Civita tensor e^,^\K, which is a genuine 4 dimen- 
sional object, and consequently with the pseudoscalar 75 
Dirac matrix. The generalization of the 75 matrix in D 
dimensions represents a problem, since it is not possi- 
ble to simultaneously retain its anticommuting and trace 
properties. In practice, the ambiguity arises when eval- 
uating a trace containing a 75 and pairs of contracted 7 
matrices and/or pairs of Dirac slashed loop momenta ^. 

To deal with a 75 matrix, several possible schemes have 
been proposed in the literature. Following the previous 
calculation of the pion transition form factor pO| , we have 
compared two of them in the present calculation. 

In the so-called naive-75 scheme u&, the anticommu- 
tation property of 75 



{7^-75} = 



(Al) 



is retained, while the cyclicity of the trace is abandoned, 
so that, for example. 



Tr [75 ^7m ^ ^ ^7'^] = (D - 
Tr [7^75 ^7m ^ ^ ^] = (2 - 



- 6) Tr [75 ^ ^ ^ 
Z?)Tr[75^^^ 



,(A2a) 
.(A2b) 



The traces obtained by cyclic permutation of the matri- 
ces 75, j^, 7^, ^, ^, ji, 7'' can be divided into two classes, 
depending on the location of 75 with respect to the con- 
tracted 7 matrices: those in which 75 is outside the con- 
tracted pair as in ( |A2a| ), and those where 7s are con- 
tracted through 75 a s in (A2b). As is seen, the result 
( A2a ) and the result ( A2b ), in which the anticommuting 
property ( Al) of 75 had to be used before the contraction 
of 7 matrices can be performed, differ by I? — 4. Con- 
sequently, if the trace is multiplied by a pole in D — 4, 
there appears a finite ambiguity in the result. 

An alternative scheme has been proposed in the origi- 
nal paper on the dimensional regularization by 't Hooft 
and Veltman ||l^, and further systematized by Breiten- 
lohner and Maison lia]. In this scheme, to which we refer 
as HV scheme, the anticommutativity of 75 is abandoned 



replaced by 








{7^^,75} = 


= for ^ = 


= o,-- 


• ,4 


[7m '75] = 


= for ^ = 


= 4, •• 


■ ,D 



(A3) 

For calculational purposes, it proves useful to introduce 
the following notation m^ : 



7m ^ 7m + 7m' 

^/i ^M ^ M ' 



where 



and 



r. 



D-A, 



The relation ( [A^ ) can then be written as 



7m 75 



-757u + 757m 



(A4a) 
(A4b) 

(A4c) 
(A5) 
(A6) 
(A7) 



This prescription for 75 violates the Ward identities 
and introduces "spurious" anomalies which violate chiral 
symmetry. To restore the Ward identities, finite coun- 
terterms should be added order by order in perturbation 
theory [^. In this scheme, the cyclicity of the trace is 
retained and the traces given in ( |A^ ) become 

Tr [75 ^7m ^ ^ ^7^] = Tr [7^75 ^7m ^ (^ A 

= {D-6) Tr [75 ^ ^ ^ ^] . (A8) 

As is seen, the result (| A8|), ob tained in the HV scheme, 
corresponds to the result ( |A2a| ) , obtained using the naive- 
75 scheme ^. 

If a trace contains an even number of 75 matrices, 
then the property 7I = 1 can be used to eliminate 75 's 
from the trace, and the Ward identities are preserved if 
the naive-75 scheme ( [Al| ) is used [|l^ (the cyclicity of 
the trace is restored and the corresponding results are 
unambiguous). On the other hand, in the HV scheme 
the "spurious" anomalies can occur owing to the non- 
anticommuting property of 75 ^. 

As for the traces containing an odd number of 75 ma- 
trices, we are left with the above mentioned ambiguities 
in the results. 



^ The presence of a pair of Dirac slashed loop momenta leads in 
fact to the appearance of a pair of contracted 7 matrices, since 



the loop integration 



'I 



I'^r 



{2n)^ {Denominator) 



Jo 



transforms I'^V'^k, ■ ■ -"t-T into 7k • • ■ 'y'^ig + ■ 



Apart from the 



contracting 7 matrices and pairs of Dirac slashed loop momenta, 
the rest of the trace elements could be treated as 4 dimensional, 
so their {anti)commutation with 75 does not make a difference. 



As stated in [|l7[, should we allow to anticommute 75 before 
continuation_ta D dimensions and after that use the HV scheme 
defined by dASl), different results would emerge and ambiguity 
would rea pp ear. We would again obtain two classes of results: 
the result (AS) for 75 outside the contracting pair of 7 matrices, 



and the result ( A2t ) for anticommuting the matrix 75 in between 
the contracting 7s. In this sense, the HV scheme would also lead 
to ambiguous results. 

Again, should we allow to anticommute 75 before continuing to 
D dimensions and then to use the HV scheme, different classes 
of results could bo obtained, depending on the position of 75 s. 
Hence, in this scheme, the ambiguity would still be present. 
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2. The 75 ambiguity in the 7*7 -^ qq amplitude 

There are two approaches one can take in order to re- 
solve the 75 ambiguity problem in practical calculations 
(such as the calculation of the hard-scattering amplitude 
of the pion transition form factor). For each diagram, one 
can determine the way of manipulating the 75 matrices 



so that the Ward identities are preserved ||9|, 10 . Alter- 
natively, one can perform the calculation using the HV 
scheme, and then introduce an additional renormaliza- 
tion constant which eliminates the "spurious" anomalies 
introduced by this prescription. In |10 , the counterterm 



at the NLO order was calculated for the special case of 
hard-scattering amplitude already convoluted with the 
asymptotic distribution amplitude. 

We have calculated the LO, NLO, and 7iy-proportional 
NNLO terms to the hard-scattering amplitude using both 



the naive and the HV prescription for 75. By combining 
these results with the results for the distribution ampli- 
tude obtained in the same order, ambiguity in the naive- 
75 scheme has been resolved and the HV renormalization 
constant determined. 

Although our end NLO result for the hard-scattering 
amplitude agrees with the result given in [nO| (determined 
up to 0{e'^)), the same is not true for the contributions 
of the individual one-loop diagrams of Fig. when cal- 
culated using the naive prescription for 75. 

Namely, owing to the fact that the relative position of 
75 with respect to the Dirac slashed loop momenta was 
ignored, classes of terms appearing in the contribution of 
individual diagrams were omitted in pfl . 

For example, the trace corresponding to the diagram 
yl23 of Fig. H leads to three different results depending 
on the position of 75: 



J 



Tr [75 pr {x r~ i) 7A {x r- i- /) i" (-(i - x) r- /) 7I 
^ Tr [7S5 /'T^ (^ r~ i) 7A {x r- i- /) i"" (-(1 - x) r- /)] 
+ Tr [(-(1 - x) />- /) 7^5 r^^ (x r- i) IX [x r- i- /) 7I 

I 



(A9a) 
(A9b) 
(A9c) 



In Rc f. ||1 [1 , however, only the results for traces (A9a) 
and (|A9b| ) were given. To conveniently describe the 
75 ambiguity present in the contributions of A23 and 
other diagrams, we have introduced the parameters 5 
and 5' . Thus, i5 = corresponds to the situation where 
75 lies o utside the contracting 7s and pairs of /Is, as 
in (A9a). Corresponding to the case when 7's are con- 



tracted through the 75 matrix is (5 = 1. F or 75 not being 
placed between the pairs of /Is, as in (A9t), 5' = , 
while for 75 placed between the pairs /Is, as in ( |A9c| ), 
S' = I. The contributions obtained using the HV sc heme 
correspond to S ^ 0. The contributions defined in (4.6a) 
of the individual NLO diagrams of Fig. Q can then be 
parameterized by 



J 



Ta11{22) — 



Ta33 — 



-^{[rjuv + 2S]-l)-^{[7jjn-2S] + l) 



-[r]uv + '2S]+ln{l-u)-l 



0(6) 



0(6) 



(AlO) 
(All) 



Ta23 — 
Ta13 = 



[rjuv + 2(5(1 + 5')]+ 2riiR + ln(l - u) - 4 



[r]uv + 26(1 + S')]+ 2r^iR 1 + - ln(l - u) 

' u 



0(6) 



(A12) 



+ - ln^(l -u) ln(l - w) - 4 

u u 



0{e) 



(A13) 



Ta12 = 



1 — V 

-2 \n{l-u)[r^iR + 26{l + 5')] 

u 

\n^{\ -u) + 10 ln(l - u) 

u u 



0(6) 



(A14) 



r 



where to facihtate the comparison with the results of Rcf. ||l0[, we have introduced rjuy — 1/e + In^^/Q^, rjiR — 
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-l/e-ln/iVQ^ and 1/e == 1/e - 7 + ln(47r). In 
results correspondin g to S' — 1 in ( |A12| ) and (^ 
well as to 5' = in (|A14|) , were omitted. 

It was argued in ||lCf| that since the quark propagator 
and the photon vertex corrections were related by the 
Ward identity of QED, they should be calculated as if 
they were not part of the trace with 75. This then de- 
termines the choice 5 = in diagrams All, A22, A33, 
A23, and A13. The remaining ambiguity associated with 
the IR (collinear) pole in diagram A12 was resolved by 
repeating the calculation in the "equal-mass regulariza- 
tion" , in which the quark and the gluon were given the 
same small mass m. It was found that the correct choice 
for diagram A12 corresponded to the result obtained by 
contracting the 7s through 75 , and for 75 placed between 
the /s (although the latter requirement was not stressed 
in [1^), that is, ior S = 1 and 6' — I. 

We have confirmed these choices by our calculation. 
In presenting our results in Sec. IV, we have already 
adopted the 6 = choice (or, equivalently, the HV re- 
sults) for diagrams All, A22, A33, A23, and A13, as 
well as for the diagrams obtained from these by insert- 
ing the quark vacuum polarization loop. The ambiguity 
present in diagram A12 and corresponding vacuum po- 
larization diagram, is parameterized only by S, with 6' 
taken as 1. 

A remark concerning the prese nta tion of the 7*7 —f qq 
amplitude in the form given by {2A) is in order. 

The multiplication of the amplitude by the factor 
l/(ie^) N Cfj^^a'fj'P" Q^ and the contraction of the Levi- 
Civita tensors before the loop integrals are evaluated 
simplifies the calculation (since then at most three-point 
Feynman integrals appear). Since the generalization of 
the Levi-Civita tensors to D dimensions is not unique, 
caution should be exercised when contracting these ten- 
sors. When using the HV scheme, the relation m^ 



,A'lA'2M3A'4 , 



7reS4 «=i 



should be employed, in which case N = 2/Q'^. Note that 

the loop integrals with I terms in the numerator of the in- 
tegrand will be encountered. Only the one-loop integrals 

containing I terms are different from the corresponding 
counterparts with / in place of I, while the presence of the 
I ■ ■ -l terms does not alter the usual results. When 



using the naive-75 scheme, the ffj,' ., from ( A15 ) can be 



replaced by g^;^^j, and then TV = 2/{{D-2){D-3)) 2/Q'^. 
We have checked our results by evaluating the contribu- 
tions with and without the Levi-Civita contraction. The 
results agree. 

APPENDIX B: FEYNMAN RULES FOR THE 

PERTURBATIVELY CALCULABLE PART OF 

THE DA 



In this sect ion we list the Feynman rules for the 0(u, t) 
operator (3.6) rederived following ]20| . 

In Section II the distribution amplitude 0(u, t) for a 
state compo sed of a free quark and antiquark has been 
introduced (3.6) 



4>{u,t) 



dz 

~2^ 

X ( 



J(u-(\~u))z- 11 



vI/(-z)I-^fi*(^) 



qq;t 



iw. 



with z+ = zj^ = 0. The quark and antiquark carry the 
momenta tP and (1 — t)P, respectively, and the frame in 
which P+ = po 4- P^ = 1, P- =P^ - P3 
has been chosen. The path-ordered factor fi ( |3.3|) 



0^j_ = 



fi 



exp < ig 



, I dsA+{ 



zs)z^ /2 



makes (/)(u, t) gauge invariant, and it can be expanded in 
perturbation theory as 



n 



n=0 



¥/:n-.-( 



Skz)z /2. 



(Bl) 



The path-ordering is immaterial since A'^ fields com- 
mute. In the light-cone gauge (A+ = 0), this operator is 
unity, but generally (for example, in the Feynman gauge 
we are using) int rod uces extra diagrams. The n-th or- 
der term in the ( |Bl[ ) series will correspond to diagrams 
with n gluon lines attached to the operator vertex. By 
inserting the term 



[1 - sk)e{i + sk) ^ I 



J 271 r^ 



(B2) 



the limits of the integration in (Bl) are changed, and the 
function 0(m, t) takes the form 
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(M) = / '-^ 



^i{u-{l-u))z- /2 



x^M:j^i(e--e- 



n=0 
X (0 



27r r 



/oo 
-OO 



m{-z)l-^A+{s,z)---A+{s^z)^{z) 



m)\^ 



/Nr, 



(B3) 



r 



The fields in (B3), along with the standard quark-gluon 
interaction insertions, are contracted in all possible ways, 
yielding different Feynman diagrams contributions. Here 
we list the Feynman rules derived from (B3). 



The operator ( 



^i- 



7"*" 75 
2^/2 



^(z) fromEq. (|B^) will 



be represented by the crossed circle (^ in Feynman dia- 
grams. The (g) vertex has a quark line entering and leav- 
ing the (S) vertex, and an arbitrary number of attached 
gluon lines (A+ field can be contracted only with A~ , so 
the gluon line does not re-enter the vertex). 

The general form of the (8> vertex with n gluons at- 
tached to it is given by 



7"^ 75 
2^2 



lr=-^--'^ric»);i(i 



(B4a) 



where 



S — u — {1 — u) — k^ — k2 — 






(B4b) 



and (7i , • • • , (?n are 4- momenta of the gluons entering the 
circle, while k\ and k^ are the 4-momenta of the quarks 
entering and leaving the circle, respectively. The form 
given in pO| is slightly different and incomplete regarding 
the gluon 4-momenta sign convention. In the special case 
when there are no gluons attached to the ® vertex, and 
the notation k\ = k and fc2 = k—P is used, the expression 
(B4) takes the form 



7'^75 
2V2 



(5(w-fc+). 



(B5) 



We will not specify here the usual Feynman rules al- 
ready used in the calculation of the hard-scattering am- 
plitude 

Let us just note that similarly to the calculation of the 
hard-scattering amplitude, the correct spin and parity 
state of the qq state has been projected by multiplying 
the diagrams by 



757" 



_ 75/' 

2^/2 V2 



(B8) 



and taking the trace. The color singlet nature of the qq 
state is taken into account by including the factor 









(B9) 



and consequently the trace over color indices must be 
taken. The result should be multiplied by an extra factor 
1/V^ (see Eq. (|B^ )), which takes into account that 
(j){u,t) is normalized to give the LO result 6{u — t) (i.e., 
normalized to 1). Finally, note that we are investigating 
the meson flavor nonsinglet distribution amplitude. 



APPENDIX C: ON THE COUPLING CONSTANT 
RENORMALIZATION 

In this section we briefly resume the relevant ingre- 
dients of the coupling constant renormalization, which 
are frequently obscured in practical applications found 
in the litera ture , and we also introduce our specific rep- 
resentation (p^) for the MS renormalization used in this 
calculation. 



The gluon propagator for the gluon attached to the (E) 
vertex (stemming from the A^ A^ contraction) takes the 
form 



9t 



q-^ -\- irj 



(B6) 



so after the g^^'^ contraction, the gluon-quark vertex for 
a gluon attached to the vertex becomes 



_A„ 



-Wl 



(B7) 



1. Coupling constant renormalization 

It is well known that in I? = 4 dimensions, the QCD 
coupling strength g is a. dimensionlcss quantity [g] = M° 
{M denotes the mass unit), while in Z? — 4 — 2e di- 
mensions [g] = M"^. Obviously, the dimension of the 
"bare"couphng constant as, related to "bare" g hy g"^ — 
4:7Tas, corresponds to [as] = M^'^. The renormalized 
coupling, i.e., the running coupling as{^J-'^), is, natu- 
rally, a dimensionless quantity and the coupling constant 
renormalization introduces the additional scale /i^ , whose 
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presence balances the dimensions in the renormahzation 
equation as M"^*^ = otsilJ^'^) Z^- The scale introduced 
by the renormahzation of the coupling constant is of- 
ten referred to as renormahzation (or coupling constant) 
scale and denoted by jjl^, while the renormahzation of 
the coupling constant in the MS scheme (the simplest 
renormahzation scheme) reads 



as = Mfl asilAi) 1 



"s(Mfl,)^ 1 



47r 



/3o- + 0(a^5) , (CI) 



where (3q = 11 — 2/3 Uf. 

In practical calculations, the additional scale /i^ is of- 
ten introduced before the coupling constant renormahza- 
tion (the presence of y? in Feynman integrals is quite 
standard), which corresponds to introducing the dimen- 
sionless {[as] — M") "bare" coupling constant as related 
to ghy g^ = Airas /i^*^, and consequently as ~ ^~^'^ . The 
"bare" and renormalized coupling constants are then re- 
lated by as = as(/x^) Za, and, if we choose to renormal- 
ize the coupling at the renormahzation scale /x^ different 
from the scale fj? introduced by regularization, the renor- 
mahzation in MS scheme reads 



as 



f^l 



^] ^sip^l) 



4^ Po- + O[as) 



(C2) 



We adopt the latter definition of the "bare" coupling con- 
stant as dimensionless quantity, and neglect the 0(a|) 
terms in further considerations. 



From (C2) one trivially obtains the scale changing re- 
lation 



as{i^^ 



mI 



^ ) -siA) 



(C3a) 
K4,) (l + ^4?^/3o In ^) + 0(e) , (C3b) 

the (3 function 



as{ 



(3{as{fi'),e) = f^^^asifi^) 



d_ 



-easifJ- ) 



2^ (^l(P^) 



in 



Po , (C4) 



and, consequently, the running coupling 

, 2< 47r 47r 

asifi 



Po €) -lU /5oln 



0(e). (C5) 



'7^ 



It is safe to ignore the 0(e) terms in (C5), since the ex- 
pression for the running coupling is usually introduced 
after all singularities have been removed (by renormal- 
ization and/or factorization). Generally, one cannot ne- 
glect the 0(e) terms in (|C2|) or (C3) (neither in order as 



nor in a|), and the use of the compact forms (C2) and 



( |C3a| ) is preferred until all singularities are removec 



2. Renormalization schemes 

Next we turn to the choice of the renormalization 
scheme and the representations of that choice. 

One can introduce different renormalization schemes 
by modifying the renormalization constant Za to Za as 
a function of 



f{e) = l + ef 



(1) 



'f 



(C6) 



where /(e) defines the renormalization scheme choice, 
while as = as(/^^)-^Q- Equation (C2) is then gener- 
alized to 



"S = ( -f ) asinn) 



An e 



(C7) 



and the choice /(e) = 1 = fMs{^) corresponds to the 
MS scheme. Equations (C2-C5) get modified by /3o -^ 

Pofie). 

The MS scheme is defined by 



/l7sW = l + £(-7 + ln(4^)) + 0(e2) 



(C8) 



and, to the order we are calculating, the 0{e^) terms are 
arbitrary. Different definitions can be found in the lit- 
erature: the original definition of the MS scheme PQ] 

/-jyg(e) = 1 + e(— 7 -I- ln(47r)) or the choices f-Msi^) = 
exp(e(-_7_+_ln47r)) and /^^(e) = (47r)7r(l-e) (for ex- 
, respectively), which mimic the 
ln"47r, (tt^)" proportional terms 



ample, |27 



and g 

e dependence of the 7' 
introduced by dimensional regularization. Although they 
are all valid choices leading to the same (MS) result 
(after renormalization and factorization of singularities), 
they unnecessarily complicate the calculation, since the 
intermediate results should be expanded over e and since 



* One of the motivations for this short summary on the coupUng 
constant renormalization was the appearance of different forms of 
ag renormalization equations in the literature. The quite often 
used form 



as = asif^R) 






i-ln4 



represents an "effective" expression, which can be strictly used 
only for calculations in which there are no singularities apart 
from those that get renormalized by this coupling constant renor- 
malization. The curious looking 05 renormalization procedure 
given in, for example, |27| , presumably represents an attempt 
to generalize the above given "effective" form to all orders in e. 
In the presence of additional UV or IR singularities, the terms 
containing both scales /x|, and ^j? remain after the 05 renormal- 
ization, which is clearly inconsistent. The final finite results are 
correct since these curious looking terms are moved to renormal- 
ization and/or factorization constants. 
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they do not contain the (F) functions originally intro- 
duced by dimensional regularization. A more appropri- 
ate choice would be the one that contains combinations 
of Fs that naturally emerge in the calculation. For this 
calculation, in which both UV and IR singularities were 
regularized by di mens ional regu larization, the appropri- 
ate choice is (see (|5.lD and (|5.4D) 



/ms W = e r(e)F(l - 6; 



r(i-e) 

F(l - 2e) 



(47r) 



eF 



^"^ (e) 
(C9) 



while, for example, for the calculation in which only UV 
singularities were regularized by dimensional regulariza- 
tion, the choice /^^^(e) = eF(e) would be appropriate. 

Alternatively, we can represent the dependence of the 
coupling constant renormalization on the renormaliza- 



tion scheme by as fie) 



') Za (see, for example, 



28[). The generalization of (|C2|) is then given by 



as 



f)\m)-' as{^.l) (l 



An e 



(CIO) 



and it rep resents a representation alternative to 
Equations (|C3|-C5) are valid for the (CIC) representation. 
Since /3o (and (3i) does not dependent on the choice of 
the renormalization scheme, the renormalization scheme 
and scale dependence can be described, to the order we 
are calculating, only by one parameter, (for example, the 
scale). The renormalization scheme and the renormaliza- 
tion scale are treated on the same footing in representa- 
tion ( |C10|) and their equivalence is explicit. For example, 
by substituting ^l = fl^ (/(e))7, A^ = A^^^ (/(e))i 



into (CIO), the renormalization equation ( |C2| ) in the MS 
scheme is obtained. 



The another advantage of the ( CIC ) representation is 
that after the coupling constant renormalization is per- 
formed, the dimensional parameter e remains as the only 
artifact of dime nsio nal regularization. Consequently, 
when using Eq. (C7), the renormalization and factoriza- 
tion constants Z^^rem Zm,coI contain 7", In" 47r, (tt^)" 
proportional terms apart from the 4r poles. In contrast, 
when using Eq. (CIC) ZM.rem Zm,coI contain only sim- 
ple 4r poles. 



In this work the representation ( CIO ) along with the 
definition (C£) of the MS scheme is used. 
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